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This paper shows how to obtain non-rigorous mathematical control over models of loosely coupled 
disordered grains; it lays groundwork for rigorous proofs and provides new information about saddle 
point structure and perturbative corrections. Both the Wegner model and a variant due to Disertori 
are transformed to matrix models which are similar to the supersymmetric model of disorder, having 
two matrices and Q'' which correspond to the two bosonic sectors of the SUSY matrix. However 
the Grassman (fermionic) sector of the SUSY matrix is omitted, and compensated by a spectral 
determinant. The transformation is exact for Disertori's model, while for the Wegner model it in- 
volves an integral which can be approximated while maintaining mathematical control. After this 
transformation a saddle point approximation is used to integrate the matrix eigenvalues, resulting in 
a sigma model. Previous derivations of sigma models of disorder assumed a spatially uniform saddle 
point independent of the Goldstone bosons and found that corrections are well controlled in the 
large A'' limit. This paper takes into account spatial fluctuations of the Goldstone bosons and finds 
that corrections to the sigma model approximation in extended systems are controlled by powers 
of the inverse conductance 1/g. The sigma model approximation is well controlled only well within 
the localized regime and far from the band edge. In other treatments energy band information is 
repackaged within phenomenological parameters; here this information is preserved explicitly, with 
clear prescriptions for incorporating it into calculations of observables. After performing the sigma 
model approximation this paper specializes to the weak localization regime, where the kinetics are 
diffusively small but at the same time dominate the disorder. In this regime Disertori's model ex- 
hibits remarkable simplifications and is completely controlled by perturbative expansions in various 
small parameters. The Wegner model likely can also be controlled but this would require further 
analysis of the saddle point equations and the determinant. The simplified weak localization model 
seems to be equivalent to the SUSY sigma model. The standard weak localization results of the 
supersymmetric sigma model, including anomalously localized states, are reproduced and extended. 

PACS numbers: 72.15.Rn, 64.60.De, 72.20.-i, 64.60.Cn 



Fifty-one years ago Anderson argued that materials with sufficiently strong disorder - i.e. complex irregularities 
at a fine scale - do not conduclii. After many theoretical developments, rigorous mathematical control has been 
obtained only in the fully localized phase-, and in the tails of the conduction band. Wegner's model of weakly coupled 
disordered grains with N orbitals per grain'^, along with the supersymmetric sigma model^^ which can be derived from 
it, are believed to capture the essential physics of both disordered conduction and localization, but have remained 
mathematically tractable only in the conducting phase with small disorder, or in special geometries. Even here full 
mathematical rigor is not preserved except in D = {0, 1} dimensions. 

Wegner's model defines an ensemble of random Hermitian'^^ Hamiltonians and prescribes that observables be aver- 
aged over that ensemble. It is conventional to adopt a field theory approach oriented toward computing this theory's 
observables using certain generating functions which can be written as path integrals. For instance, the energy level 
correlator i?2(£'i, E2) may be written as the second derivative with respect to i? of a generating function: 



The ipvni vector contains 2N complex Grassman variables at each site, while Syni contains 2A'' complex scalars at 
each site, n specifies the orbital, v specifies the lattice site, and the two possible values of i = {1,2} correspond to 
the two energy levels, whose energies are given by -E^. L is diagonal sign matrix introduced for convergence reasons. 

When E is proportional to the identity, the Wegner model has an exact symmetry under unitary transformations 
tpi — » Uii^^ir^, and another exact matrix symmetry under hyperbolic transformations of S. Therefore one expects 
that if the energy level splitting uj — Ei — E2 is small enough, then the physical degrees of freedom are not 2A^- 
component vectors but instead two 2x2 matrices: one Hermitian and the other a member of the pseudo-Hermitian 
group corresponding to hyperbolic transformations^. 



I. INTRODUCTION 




2 



Efetov pursued a program of first averaging over all Hamiltonians and then making an exact change of variables 
from the original vectors to matrices. His result has at each site a single 4x4 graded matrix Q which contains the 
two required Hermitian and pseudo-Hermitian matrices, along with additional Grassman variables in the off-diagonal 
blocks^. A non-local and nonlinear action ^■Tr{Q^) + NTr{QE) + NTr{ln{Q + ek)) controls the graded matrices^. 
The N multiplying the action invites a saddle point approximation which would constrain Q's eigenvalues and leave 
only angular variables as dynamical degrees of freedom. Although a correct application of the saddle point technique 
would require taking into account fluctuations in the saddle point induced by the fluctuations in the angular variables, 
until now all sigma model derivations have neglected such fluctuations and have assumed that the saddle point is 
spatially uniform^i^i^ii. It has also been conventional during this procedure to avoid explicit consideration of the 
details of the band structure, and to instead introduce phenomenological constants. After fixing the eigenvalues 
and making a Taylor series expansion of the logarithm in powers of Q's fluctuations, Efetov arrived at the famous 
supersymmetric sigma model. 

In this article I follow the slightly different path of Fyodorov, as developed by Disertori. In the case of a single grain, 
Fyodorov showed how to make an exact transformation^ from the original Wegner model to 2 x 2 matrices and 

. Disertori followed Fyodorov very closely but considered Hamiltonians which are entirely random, and showed that 
Fyodorov's exact transformation succeeds in any dimension on a very wide class of lattices^. The principal difference 
between the two models is that Wegner's density of states is determined by the kinetic operator fc, while Disertori's 
density of states is just the semicircular distribution of random matrix theory. We analyze Disertori's model and 
consider also the Wegner model. The latter model contains a third matrix degree of freedom in addition to the 
and matrices. These W matrices exhibit spontaneous symmetry breaking and therefore can be integrated while 
maintaining mathematical control. 

The matrix models produced by Fyodorov's transformations are similar to Efetov's model: the Lagrangian's terms 
have direct correspondences to Efetov's action ^Tr{Q^) + NTr{QE) + NTr{\n{Q + ek)). However there are no 
Grassman variables in Fyodorov's matrices, and in compensation the path integral contains a determinant coupling 
with Q^. The overall effect is the same as if one had started with Efetov's model prior to any approximations 
and then somehow managed to integrate exactly the Grassman variables. This route has been until now impassable 
because Efetov's logarithm contains all powers of the Grassman variables. Fyodorov's alternative procedure bypasses 
this difficulty and is still exact in Disertori's case and well controlled in Wegner's case. 

Having converted to matrix variables, we analyze Disertori's model in detail, taking care to maintain mathematical 
control at each step. There are two principal challenges. The first is to integrate the matrix eigenvalues in a controlled 
manner, in contrast to the conventional SUSY procedure which presumes a spatially uniform saddle point. In contrast 
to previous treatments we find that the sigma model approximation to both the Wegner model and the Disertori model 
is corrected by powers of l/g not 1/N, and fails near the band edge and in the localized regime. The second challenge 
is to control the — coupling. This is a very intricate object, but simplifies drastically in the weak localization 
regime, where fluctuations in and are small and the kinetics dominate the disorder. In this regime the — 
coupling simplifies to a zero-momentum coupling, a normalization constant, and corrections proportional to 1/g. In 
the end we reproduce and extend the standard SUSY results about the two level correlator and anomalously localized 
states, including the correct normalization constant. Our results go beyond the SUSY results by incorporating explicit 
band information. Additionally we find in the two level correlator an extra term controlled by on site elements of the 
Green's function, but this term may have been simply overlooked by previous authors, or could possibly be cancelled 
by higher orders of 1/g perturbation theory. A partial analysis of the Wegner model suggests that it also agrees 
with SUSY results for observables in the weak localization regime, but completion of this analysis waits on further 
perturbation theory results about the dW integral and more precise solutions of the saddle point equations. 

Our results suggest that the simplified weak localization version of the models treated here is equivalent to the 
SUSY sigma model. Corrections to the sigma model approximation and explicit band structure information could 
likely be obtained within the SUSY formalism. What seems more doubtful is whether outside of the weak localization 
regime SUSY could correctly reproduce the present article's derivation of the correct — Q'' coupling. 

Section |TT] of this article performs exact derivations of matrix models equivalent to Wegner's and Disertori's models. 
The rest of this paper concentrates on Disertori's model, but includes asides about the Wegner model. Section IIIII 
performs the saddle point approximation, section IIVI obtains control of the — coupling in the weak localiza- 
tion regime, section |V] integrates fiuctuations in and Q^, section IVTl shows that the standard SUSY results for 
weak localization are reproduced and extended, and section IVIII reviews how we maintained mathematical control 
throughout. Lastly we briefly discuss new vistas for understanding disorder outside of the weak localization regime. 
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II. DERIVATION OF THE MATRIX MODELS 



In the next pages we will develop Fyodorov's approach^ to disordered systems. Beginning with a general model 
of loosely coupled grains, we will move by exact steps to a path integral with vector degrees of freedom and then 
to matrix degrees of freedom. The Wegner model and Disertori's model will be considered together through most of 
the derivation. Both models prescribe random ensembles of Hermitian Hamiltonian matrices modelling conduction 
through a system. The system is modelled as a graph with V sites, and we will avoid taking the continuum limit in 
order to preserve flexibility and mathematical rigor. There are N basis elements at each site, so the total basis size 
is N xV. The lower case letters n,v will denote respectively the orbital index and the position. 



A. The Observables 



The field theory approach to disordered systems computes generating functions which are tailor made for specific 
observables; therefore one must start by choosing the observable to be computed. In this paper we will calculate the 
averaged two point correlator, averaging over the ensemble of Hamiltonians: 



Tr{5{Ei^ H))Tr{5{E2- H)) 



R2iE,,E2)^ ' ' ^, p{E)^TridiE-H)) (1) 

p{Ei) p{E2) 

S in this context is the matrix delta function. There are also spatially resolved versions of these quantities; the first 
line of the next equation writes them in terms of matrix delta functions, while the second line gives them in terms of 
H's eigenvalues e and eigenfunctions ipe- 



P(E,V) = y^{5{E - H))nnvv, R2{Ei,E2,Vi,V2) = 

P[El,Vi) p[E2,V2) 



-jjr-, Sr.tj, W/IM2 T,fj, r ,_ J:....HEi-e,)diE2-e2)\{AAviW\{AMW 
P[E,v) = 2^6{E - e)\{i^4v)\ , R2{Ei,E2,vi,V2) ^ — - (2) 

^ p{Ei,Vl) p[E2,V2) 

In the zero-dimensional gaussian unitary ensemble R2 has been calculated using many different techniques. Taking 
e as the characteristic scale of the disorder, it is 



TT^uj^{p{E)Y Tre V 4e^ 



R2{E + ^/2, E - C./2) = 5{u:p{E)) + 1 - ^ \,-^': , p{E) = -\/ 1 - 7^ (3) 



B. Green's Functions and their Generating Functions 

In the field theory approach one rewrites the delta functions in terms of Green's functions. It is important to 
distinguish between the advanced and retarded Green's functions: the former is Ga[E) = {E — iv — H)^^ , while 
the latter is Gr{E) = {E -\- iv — H) . v is infinitesimally small, positive, and destined to be set to zero sometime 
during the calculation. The scalar 5 function can be represented as ^ limi,^o Im{{E — iv — E) ); therefore p{E) = 
i lim,y^o I'm{Tr{GA{E))). Similarly, the two point correlator can be extracted from product of two Green's functions: 

p{E)p{E) = \\iinIm{Tr{GA{E)))Im{Tr{GA{E))) 

= -^\im[iTr{GA{E))-Tr{GR{E))) x {Tr{GA{E)) ~ Tr{GR{E)))] 

= lim[Tr{GA{E))Tr{GA{E)) + Tr{GR{E))Tr{GB.{E)) 

47r i^^o 

~Tr{GA{E))Tr{GR{E)) ~ Tr{GR{E))Tr{GA{E))] (4) 
= :^^^5^=^£^, Rj,A = \\u,Tr{GR{E))Tr{GA{E)), Raa = \imTr{GA{E))Tr{GA{E)) 



In general, calculating an n-point correlator requires averaging products of n Green's functions. 
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It is mathematically convenient to calculate not Green's functions but instead their generating functions. First 
consider a function which is sufficient for generating a single Green's function: 

f , _ det{Ef-H-jf) 

J is a source matrix J = Juiua^nma- has a small imaginary part v which one chooses to be either negative 

if one wants an advanced Green's function or instead positive if one wants a retarded Green's function. All other 
quantities are Hermitian. When the source J is set to zero and the two energies are set to be equal, the generating 
function Z{E^ = E^, ji = = 0) is equal to one. Using the identity det{A) = e^''^'"^'^^) one can easily prove that 

^ = Z X E„ {nv2\{E^ -H- J^)-'\nv,). 

V1V2 

One can obtain a Green's function from Z by first taking its derivative with respect to and then setting 
E-f — E^, — = 0. Similarly, one can take a derivative with respect to set E^ = E^, = = 0, multiply 
by —1, and obtain the exact same result. The identity of the derivative with the J* derivative is a Ward identity 
for this theory, a manifestation of the symmetry between the determinant in the numerator and the determinant in 
the denominator. The numerator-denominator symmetry is part of the graded matrix supersymmctry of the SUSY 
sigma model, and is very important for understanding disordered systems. 

If one wants globally averaged quantities one may take derivatives with respect to E; J is useful only for spatially 
resolved quantities. We will put a tilde over J as a reminder of this. 

The generating function for a product of I Green's functions is a product of I ratios: 

^^ nlidct(ii;/-g-j/) 

n;=iciet(4-//-j}) 

Instead of keeping track of all 21 determinants separately, it is convenient to merge them into one determinant in 
the numerator and one in the denominator. For instance, consider two matrices A and B which each inhabit a basis 
with M basis elements. One can invent a new basis with 2M basis elements - M elements for A and M for B. In 
the new basis AB = BA = and therefore det {A) x det {B) = det {A + B). This is a sleight of hand: one avoids the 
complication of having two determinants at the expense of using a basis which is twice as large. We apply this trick to 
Z ^ where each determinant inhabits a basis with NV basis elements. In order to group the numerator's determinants 
together, we create a new basis with NVI basis elements, denoted by the indices nvi. The denominator's basis is 
similar, denoted by the indices nvj. The new generating function is 

^ ^ t(£^^ H t/^) ' — ^viV2^iii2^nin2i ^ — ^j^viV2^jij2^nin2i 

= JiviV2^iii2^nin2J ^ = Jji,^i,2^jij2^nin2i = Hy^^^5i^i25n^n2i ^ = Hyj^y^Sj^j^Sn^n^ (7) 

We will consider a more general generating function with determinants in the numerator and determinants 
in the denominator: 

Z ^ n[Iidct(i;/-g- j/) ^ det(i'/ -H-P) 
rijli dct(£'| -H- Jf) dci{E^> -H-.P) 

The extra generality can be useful for obtaining a toy model of unquenched QCD, for checking the theory's constants 
and scaling, and for investigating various special cases. 

The next step is to rewrite the generating function as a path integral, converting the determinant in the numerator 
into an integral over Grassman variables and the determinant in the denominator into an integral over bosonic 
variables. If A is an Hermitian M x M matrix, tp, tp are vectors of M Grassman variables, and 5 is a vector of M 
complex bosonic variables, then 

r - - r _ ^ _ 

det{A) = {-2i)^ / d^d^pe^'l"^'!' , det'^A) = (27ri)"*^ / dS^dS^e^^^^' , d^^diP ee JJdVm^m (9) 

However the bosonic integral diverges if the imaginary part of A is not positive definite. This seemingly easy point 
causes some of the interesting mathematical intricacies in the supersymmetric theory. Using these identities, we 
rewrite the generating function: 

Z = {27T)-'''''\^/2)-'''''\^)-''^'\detLf'' J #d^d5«d5^e^^^(^''-^-^")^*e**(^'-^-^'')^ 

^n\n,2V\V23i32 — SW^{-^^i-Ej^^^nin2^viV2^jij2 (-'-0) 



5 



The Tpvni vector contains N complex Grassman variables at each site, while Synj contains PN complex scalars 
at each site. Both vectors have dimensions of [Energy]'^/"^ . L is a diagonal sign matrix introduced to ensure the 
convergence of the bosonic integrals. Its signs correspond to the choice of retarded vs. advanced Green's functions, 
and are +1 for the retarded case and —1 for the advanced case. 

One could also introduce a sign matrix similar to L in the fermionic integrals, and their guaranteed convergence 
allows one to choose any combination of signs one likes. Verbaarschot et a& explored this freedom in the context of 
the supersymmetric sigma model, and discovered that in that context one must choose the signs to be all the same, 
and thus obtain a compact representation for the fermionic variables. In the approach being developed here, one can 
choose any signs one likes, but soon one integrates out the fermions and then the fermionic sign matrix factors out 
entirely; one is again forced to use a compact representation. 



C. The Models 



We will analyze models of disordered grains with V grains and N orbitals in each grain. The Hermitian Hamiltonian 
H = eoi?o + has a deterministic part Hq and a random part H which is fully described by the second moment 
Hmn2ViV2Hn3rnvav4, = N'''-{1 - fc)„j t,^ 5„j„^ J^^ng (5« 5^2i,3 . The system geometry, including the number of dimensions 
and all other structural details, is encoded in the positive indefinite kinetic operator k and in Hq. We work in the 
diffusive limit (s|A;|s) <C 1 Vs, where |s) are the momentum basis functions. We also require that k and Hq be 
Laplacians, meaning that k\0) = and Ho\0) — 0. eq and e are the energy scales of the two operators. 

The traditional Wegner model^ is obtained by setting Hq — — and fc = 0, so that the kinetic term is deterministic 
and the random potential is local. Disertor^ proposed a more tractable model where Hq = and both the kinetics 
and the potential are random. We will analyze both models. 



D. Averaging over the Disorder 

We want to calculate the average of products of Green's functions, so the next step is to average the generating 
function Z over realizations of the random potential. One expands formula [TO] as a power series in the random 
potential H, counts pairings of H, and then substitutes the second moment for each pairing. The term to be 



averaged is exp {^i{SLHS* + i^Htp)) = J^Zo IF (^) {SLHS* + i^Hip) . Odd moments of H average to 0, while 



4| \' 

even moments average to constants times powers of the second moment 

g \(SLHS*+ijH^f] ' = e-2"'^"(^^^^^^+^, Np{2a) = (27r)-^/' / daa^'^e'i"' = {2a - 1) 

i^Q y )■ I J J 



[11) 



Inserting the formula for ff's second moment completes the averaging process; 



C = iNr^Tr.,,{SL{E''-eQHQ-j''))+'-^P{Ef^eQHQ-P)^ 

- y ~ k)y^^^Tr{tpy^'ipy^) + ^ (l - k) .y ^ y^T T {S y , y , L S y^y^ L) + 2X) 

V-LV2 VIV2 

g2 _ 

^ = y^ ^ ] (1 ~ ^)viV2'4'niViiifpn2V2ii^n2V2jl^niViji-^jl 

viV2iijinin2 

1pViV2ili2 = y^y^^ V'thJIU V'"fl»2'^flf2 1 ^VlV2jlj2 = '^T^'^^, ^nvi jl ^riV2j2 (12) 

n n 

The minus sign in the quartic ijj term was caused by the anticommutation of fermions and is important for the theory's 
convergence. The new dimensionless matrices S and ip give an early hint at this theory's matrix structure. 
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E. Hubbard-Stratonovich Conversion of the Fermionic Variables 

Up to this point we have kept company with Efetov in his development of the supersymmetric sigma model^. Now 
we part from him, taking a parallel path. Both paths convert from vector variables S and tp to matrices and then 
perform various approximations. The distance between the paths is caused by a different choice of matrix variables: 
Efetov chose graded matrices, while we follow Fyodorov in choosing ordinary matrices. 

First we do an exact Hubbard-Stratonovich transformation of the fermionic vector ip into a bosonic x Hermitian 
matrix . If 1 — fc is positive definite. 



2 

VI ,V2 



exp(f kU^^^Trii^M) = (iV/2^)^^''/'2^^^'(^^-i)/2(det(l - k))''^^' 

/N 
dQexp(--^(l - k)-XTr{Q.,,Q.,,)±NY,Tr{QfJ,)), 

ViV2 V 

dQ - ([ldQu)illdQtdQiJ ^A^y^Mix)dUdx, Q^UxU^ (13) 



UxW gives Q's decomposition into eigenvalues x and a unitary matrix U. AvdM{x) — YiiiKi^ (^n "^'2) the Van der 
Monde determinant. We apply this transformation at each site in the system's volume, with a different Q-^ = Ux^U^ 
at each point. 

Z ^ -f f dU dx (fpdipdS^dS^ Y[a1^m{x{,) exp(£) 

- Y E (1 - k)Zl.,Tr{Qf ,^Qf ± NY,Tr{QfJ,) - (1 - kU,,Tr{S,,,,LS,,,,L) + 2X) 

V\V2 V V1V2 

7 = (2^)-^^^\*/2)-^^^'W-^^^\deti)^^(iV/27r)^''^/^2^'(^'-i)^/2(det(l-fc))^''/^ (14) 

The new Lagrangian is linear in -0 and ip, so we use equation [S] to integrate these variables. At the same time we 
rescale S {2N/e)~'^/'^S to remove energy units from the integral, e = eo/e is the strength of Hq relative to the 
disorder. 

Z = dU dx^dS^dS' Ay^jv^(x{) exp(£) det(Ao(5„^„^-e(2jV)~^ (1 - fc)^i^^5'ni„iji-^ji5'*,^^,-/»iij. 



C = iNTr,.j{SL{r^& - eHa - r^j'')) - —J2 (1 ^ k),,,,Tr{S,,LSy,L) 

V1V2 



VlV2 



7 - (2^)-^^^\z)-^^^'-^^^'(detL)^^(iV/27r)^''^/^2^'(^'-i)^/2(det(l - k)f^l\2NY^'' (15) 

Fyodorov showed how to rewrite the determinant to depend on the S matrix rather than the S vector^'*. The 
following lines generalize his proofs to the Wegner model, where i/o and are nonlocal. The 5 vectors may be 
understood as an NV x l^V matrix (diagonal in w), and X]j "^'Hfij-^jXl ~ ^)fiW2*5'j^2i'2j rewritten as the 

multiple of matrices: SkS\ where kviV23i32 = (1 ^ ^)fif2'^iii2-^ji ■ Furthermore 5 may be written in its Singu- 
lar Value Decomposition Symj — WsC, where W is NV x NV and unitary, C is I^V x I^V and unitary, and 
the diagonal I^V x I^V matrix s is composed of S's singular valuesifl. With this notation the determinant is 
det(^°5„,„, - e{2N)-^S,,,,S{l - k)LS^) = (det A^)^ det(l - i{2N)-HA°)-^WsCkCh^W^). We would like to use 
the cyclic properties of the determinant to move C^s^W^^ around to the left side, but this is not allowed because C 
and s do not have the same rank as the argument of the determinant. Therefore we define augmented NV x NV 
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matrices S, and L, which are padded with ones on the diagonal and zeros everywhere else. We also define a 
diagonal matrix 9 which has ones in the first I^V diagonal entries and zeros everywhere else. One may verify that 
sC(l — fc)LCtst — sC{\ — k)L0G'^s^. This leaves us free to continue rearranging the determinant: 

(det A° f dct(l - i{2N)-^{A°)-^W~sC{l - k)Le&s'^W'^) 
= (det A^f det(l - l{2N)-^&s^W\A^)-^WsC{\ - k)LB) 
= (det A^ f det(l - 1{2N)-^S\A^)-^S{1 - k)L) 
= (det A")'^dct(l - {A\Xv,v,S.,.,n32Ln{^ - k),,v,) 

= {det A^f-'" dei{A^5,,,, - A^), A' ^ <nvovM'')iXv,vJ-^^.njMl ~ k^.v, (16) 

ilViV2 

The 9 in the second line prevents the P x xV sector of the determinant's argument from coupling with the rest 
of the basis. Between the second and third lines the operator inside the determinant moves from living in a basis of 
size X N X V basis to living in a basis of size P x x V . 



F. Conversion of the Bosonic Variables 



In the previous section we integrated out the fermionic vectors ^ and ip associated with the determinant in the 
numerator of the generating function Z . The result was a new matrix Q-^, which I will call the fermionic matrix even 
though it contains no Grassman variables. Now we will extract an x Hermitian matrix from the S vectors. 
This extraction process was done first by David, Duplantier, and Guitter— . Fyodorov re-derived their result and 
applied it in the context of random matrix theory^, and Spencer and Zirnbauer— proved the special case = N . 

We start with the singular value decomposition S = WsC and assume that N > . In contrast with 
the previous SVD, now each of the matrices S't,,W^,s^, and Cy are local. The integration measure becomes^ 
dS^dS^ = ds'^dWdCl\l^M{s'^)(dets'^)^-^\ We define the positive indefinite matrix = {2N)~^CW^C; dS^dS' = 
dQ^ dW 9 {Q'') {det Q'')^~^ . Fyodorov^ computed / dW. Adjusting for factors of two caused by Fyodorov's use of the 

integration measure dSdS* versus our use of dS^dS^ , this constant is / dW = (Hil^v-t '0 7r~^''(^''~^)/^+^''^. 

As a consequence of the original theory's use of L to make the dS integrals converge, always occurs in combination 
with L. FyodoroA*^ showed that Q^L factors into Q^L = Tx^T^^, where is diagonal and constrained by x''L > 0, 
T is a member of the pseudo- unitary hyperbolic group U (n+, and n+ and n_ are the numbers of plus and minus 
signs in L. In the special case where all of L's entries have the same sign, T is an ordinary unitary matrix. With 
the exception of this special case, U{n^,n-) is not compact, and as a consequence has parameterizations in which 
one or more of its parameters is unbounded. Fyodorov's integration measure is dQ^ = dx''dTAlr^j^.j{x^). In these 
coordinates the path integral is 



Z = 7 / dxf dUdx" dT Yl AI^m{4) ^IdMiO 

Jx<'L>Q 



e 



= / dWcicp {~iNTr„j{SL{eHo + r^J''))) dct{A°6j,j, - A^{S)) 



C = <k + - E + E Tt{T,<T~^E^) + ^Nr' E TT{U,xiulEf) 

vk vj V V 

— E ^^/^^ ~ k))viV2Tr{Uy^xl_^UljJy^xl^Ul^) + "^E ^'^i'"2^''(^i'i2^!^i^i;i^^'"2^!^2-^i'2^) 

ViV2 ViV2 

+ {N - I^)Tr{\nA°) 

7 = (2vr)-^^^\^)-^^^'-^~^'(detL)^^''(iV/27r)^''^/^2^'■(^'-l)^/2(det(l - fc))^''/^(2A^)^^^^ 

^ = ^iont>ofl(^ )iii2ViV2^VrV2jlj2-^jl{^ ~ k)v2V3, 

il ,V2 

S.,y,L = {2N)-'Cls.,WlW,,.Sy,C,,L, {2N)-^ClslC^L ^T,xIt~^ (17) 



8 



G. Disertori's Model 

Disertori^ applied Fyodorov's ideas to a model in which the Hamiltonian is entirely random; Hq = 0. She also 
required implicitly that the source J be diagonal in the position index u, which restricted the theory to calculating 
on-site elements of the Green's function G„„. With these restrictions in place, A° is local and det A° factorizes site 
by site. simplifies to Ql^Lj^{l — fc)„i„2'^iii2 j so that there is no dependence on W; the remaining integral in e^^ff 
is just (/ dW)^. 

Shifting ^Qf - li-^E^, we obtain Disertori's model: 

Z ^ if dUdx^ dTdx' Xl^lamiO ^Idui^) det(;7,,<f/„\<5,,,,<5,,,, - T,,xlj-\l - kU,J,,,,) 

Jx''L>0 \ 

C = <k + {N~ /") ^ Inx.fc + ^ TriT^J'J-^E' ?)) + zTVr ^ ^ Tr{U,xiulEf) 

vk vk V V 

— — ^ {k/{l — k))viV2Tr{Uvj^xl^Ul^Uv2xl^Ul^) + ^viV2Tr{Ty^x''^_^T^_^T^^x\,^T^^) 
V1V2 V1V2 

N-l ~V 

7 = N^''^+''''^/\AeiLf'\ W l\) 2-^'^/27r-^'(^'-i)^/2-^'^'^/2(det(l-fc))^''/' 

l=N-I'' 

^ ^^TriEfEf)^-VNl''-VNlf ^^^-^ 

This completes the sequence of exact steps which converts Disertori's model to matrix coordinates. The action 
is remarkably like the exact SUSY action ^;Tr{Q'^) + NTr{QE) + N Tr{\ii{Q + efc)); the primary difference is that 
the kinetics are displaced from the logarithm to the quadratic terms, in keeping with the model's semicircular band 
structure. 

After the shift — > — ie~^E-f one must remember that depends on E^ when making use of the fermionic 
sources. Neglect of this dependence will result in incorrect prefactors and will break the Ward identities associated 
with the numerator-denominator symmetry, even when calculating the density of states. To avoid these intricacies, 
from now on we will use prefactors instead of sources. Derivatives with respect to the local sources J produce the 
following prefactors: 

dJ^ J 
d^ d d 



dJl.djL. djL,, djf 



(N _ . Fi . . + ?-2 V . . . p2 

y±y J J'. "I'lf2-' i;ilil2 '"11211 ' / y fl l'2'l '2 71 72 f 2 fl »2»1 72 71 



Ulll"'"D2»2 """Uili ""'i;2l2 Jlj2 

F^ = {Ux^U^f\ F^ = {Ux^U^ -Tx''T-\l-k)f' 

^ -iNr\Tx'T-%,,, -^r^, (19) 

J t6 ^ dJ^ dp dJ'^dJ^ 

When one is interested in calculating off-site elements of the Green's function G^^^^ these prefactors require modi- 
fication to adjust for non-local sources. If one uses bosonic sources, must have an off-diagonal element connecting 
vi and V2. Therefore the integral in >Ce// does not decouple from the Wy^ integral. The correction to the path 
integral looks like 



( j dWy J dWidW2eM-^^''NTr,={,^^,^^{SLj'')), 

S,,,,L = {2N)-^Cls,,WlW,,s,,Cy,L, {2Ny^ClslC,L = T.xIt-^ (20) 
Evaluation of this integral lies outside the scope of this paper. 



H. The Wegner Model 



In the Wegner model the random potential is purely local and the Laplacian is deterministic. Within this paper's 
formalism it can be obtained by setting fc = and Hq equal to the Laplacian. I take the liberty of renaming Hq = k. 
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After absorbing into Q-f , 

dxf dUdx" dT dW n ^ldM{4) KdMiO detiA%,,,, - A\S)) 



x>'L>0 



C = -^^xlf,+iNr^^Tr{U^xlUlEf) + {N - I^)Tr^,{\n{Uxfu'^ -lek-ir^jf)) 



2 

vk 



iNr^ Tr{T,xlT-^&) - iNTr,^ {SL{ek + r^j'')) + (iV - /'') ^ In 



x''- 



7 = (2^)-^^^\z)-^^^'-^^^'(detL)^^\7V/2^)^''^/^2^''(^'-i)^/2(det(l - k)f^ /\2NY«'\^^'-(^' ^'^ 

-^ViV2iii2 — Qiiixi2^'"iV2 ~ ^^^iii2^viV2 ~ "^^iii^^ J^i ^ = ^ioiiUoUi ) iii2ViV2^'"i^2jlj2-^jli 

SviV2-^ — ^Vi'^Vi^^^^^^V2^V2^V2^^ Cy'^yCyL — TyX^Ty (21) 

This action is similar to the exact SUSY action ^Tr{Q'^) + NTr{QE) + NTr{ln{Q + efc)), where Q contains an 
implicit I. However Q^'s kinetics are hidden in the dW integral. This integral is similar to a logarithm: its derivatives 
contain Green's functions, and its contribution to the Lagrangian is a logarithm at leading order. Moreover the Ward 
identity '^^(^'=f'^=^) = Q ensures that the dW integral has the same physics as the Q-^ logarithm in equation [211 
The Q-^ logarithm corresponds to multiplying the path integral by {Q-^ — iek)^~^ ; when ek ^ this multiplier 

scales with e^~^ , and in the opposite case ek <C Q-^ it is constant with respect to e. The Ward identity ensures that 
this multiplier must be cancelled by some other contribution, and the only term which can fill this role is the dW 
integral. Therefore the dW integral must scale as ~^ when ek 3> Q-^ and be roughly constant with respect to e 
when ek . 

The Wegner model and Disertori's model are very much alike. Section IIIIDI examines the similarities of the two 
Lagrangians and their saddle point approximations, which work out roughly the same. The principal difference is 
that Disertori's energy band is the semicircular one of random matrix theory, while Wegner's energy band can have 
any shape at all, according to the kinetic operator which one chooses. Another superficial difference is that Wegner's 
model allows easy tuning of the relative strength of the kinetic operator, while in Disertori's model a single energy 
scale controls both the kinetics and the disorder. This is merely superficial: N^/'^k can be used to tune Disertori's 
model. Section flV Dl discusses the — determinant, which is wholely responsible for interactions between Qf 
and Q'^. In both models the argument of the determinant is the sum of the kinetic operator and Q-^^ , minus a term 
proportional to Q''L; it is Q^Lk + — Q^L in Disertori's case and —lek + Qf — A^{{A^)'^^ ■ SL) in Wegner's case. 
This form describes competition of the kinetics versus and Q^. We will concentrate on the case where k dominates 
the determinant, otherwise known as the weak localization regime. In this regime at leading order Wegner's — Q'' 
coupling is the same as Disertori's coupling, and both are independent of W . 

We turn to the dW integration. The Wegner model depends only on the first columns of W . Taken in isolation 
from the rest of W ^ these columns form the Stiefel manifold V/6(C^), which is defined as the set of all possible 
combinations of orthonormal complex A'^-vectors. This manifold forms an homogeneous (continuous and transitive) 
space corresponding to the unitary group; Vib{G'^) = U{N)/U{N — I^). Integration over the other N — columns 
in W produces a constant. 

The dW integral exhibits an exact continuous symmetry under global rotations of W , and is governed by a kinetic 
term —iNeTryj{SLk) which regulates fluctuations in W . Therefore Ws symmetry is spontaneously broken in I? > 2 
dimensions, and in small volumes displays an effective SSB even in D = {1,2} dimensions''^. In the SSB phase W takes 
on a spatially uniform value, with small fluctuations. Therefore we will do perturbation theory in W^s fluctuations. 

We can assess the validity of the W perturbation theory, and of the SSB assumption itself, by calculating corrections 
to the leading order results. As an example, consider the P = 1 case where the Stiefel manifold is a sphere and W 

is a unit vector e. The free energy F is proportional to In / c^e exp(— leA^ ^^^^^ ky2Vi\J'x^ ej^ • e^j^y^x^). At leading 
order F scales with the number of degrees of freedom, which is 2A^ — 1. The first correction Fi is proportional to 
^1)11)2^^11)2' where G is the bare propagator ((ej^ • Ji)(j2 • 61)2))- The numerical value of this correction 
may be calculated either numerically or by analytic approximations. When SSB occurs G — {Nex^k)^_^y^6j-^j2, and 
the relative strength of the first correction to the free energy is controlled by {Nx^ek)~^ . We will find that the saddle 
point value of x^ is independent of N. In the N ^ 00 limit SSB can be frustrated only when D — {1,2} and only 
when V ^ 00, in which case the propagator's diagonal elements diverge^^. 

In the SSB phase prefactors are not affected at leading order by the W fiuctuations, so we move the bosonic source 
J* outside of the dW integral. This step is exactly correct if one calculates only local observables like Gw If one 
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wants non-local results then perturbative corrections may be calculated in powers of the bare propagator. In the weak 
localization regime the — determinant is almost entirely independent of S and is best treated as a prefactor, 
so we move it also outside of the dW integral. Here too perturbative corrections may be calculated if desired. 

We assume SSB and calculate the dW integral at leading order. We parameterize as a unitary matrix; i.e. 
W = , where Y — —V^ is anti-Hermitian. If one calculated corrections one would find an additonal 0(S) term 
added to £i's h — h, where S cx iV is Y's self-energy. 

Yl (xlj)'^^^'' J dWexp {-iNTr^j{SL{ek + r^P))) det{A°Sj,j, - A\S)) 

vj 

= {2T:/iNf'''^'^-^''/^'^dct{A°Sj,j, - A\S)) e^'+^\ Ci = -{N - /2)Tr{\n{h - h)) + [N - I^) ^ In 
£2 = -iNTr,j{SL{ek + r\j'') +^-Y,hl-{h- h)-^ ■ K) 

V3 Vi 

Sv.v^L ^ Clsv^Sv^Cv^L, ClslCvL = nxlT-\ (Y^Y) = {1 - l' /{2N))iie)-\h - h + tJ:)-' (22) 
This simplifies tremendously when s and C are spatially uniform and efc 3> 1: 

£1-1- £2 = -tNr^TrvjiSLj'')- {N ~ l''/2)Tr\n{Lek)^ l''/2Y,^nxlj (23) 



VJ 

I 



VJ 



When efc <C 1 the correct result for Ci + £2 retains the original factor of {N — ^^^^vj^'^^^j- "^^^ W ~ 
parameterization can not be used for other computations because it mixes the Stiefel manifold with the rest of U (N) . 
E, (Wtpy^^ and corrections to the dW integral require a parameterization adapted specifically to small fluctuations 
on the Stiefel manifold. In the case of = 1 the e unit vector parameterization produces easy results for the bare 
propagator and the leading order self-energy; G oc {eNx^k)~^ , E(s) cx J2si^o k~^{s — si) fc(si). Yet even here (VF'I'W) 
contains all even powers of e's components. Further consideration of these issues lies outside the scope of this paper. 
We will focus mainly on Disertori's model, and will at intervals discuss how how similar results can be obtained from 
the Wegner model. 



I. Is the — Q'' Model Supersymmetric? 



It is important to distinguish between different meanings of the term supersymmetry: 

• Both the Wegner model and Disertori's model possess continuous symmetries associated with global rotations 
in Qf UqQ^Uq and Q''L TqQ'^LTq^ , SL TqSLTq^; these symmetries are exact when E is proportional 
to the identity. 

• We have already discussed another symmetry which is always exact: the symmetry between the determinants 
in the numerator and the determinants in the denominator. The numerator-denominator symmetry is not a 
continuous symmetry. Instead it manifests itself in any number of Ward identities, including ones relating Q** 
observables to observables. 

• In the supersymmetric sigma model there is a single degree of freedom, the graded matrix Q which contains both 
Grassman and scalar variables. The scalar sector is roughly equivalent to this paper's and matrices. In this 
model both the rotational and the numerator-denominator symmetries are subsumed into a graded continuous 
symmetry connected to Q's rotations in a graded group; this is the traditional meaning of supersymmetry. The 
supersymmetric model's numerous successes, including its noteworthy non-perturbative results, are attributed 
to the graded symmetry. 

Fyodorov's conversion to — coordinates has been discounted because it does not maintain a unified approach; 
in particular there is no symmetry transforming variables into variables or vice versa. Despite the exactitude of 
its derivation, one fears that without the graded symmetry's protection later approximations will produce uncontrolled 
results. The obvious differences between and Q*" make this fear more tangible. In Disertori's case the Lagrangian 
has stray 1 — fc factors in the — determinant and in Q-'^'s kinetics. The observable (equation [19]) is written 
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in terms of both {Q-^)~^ and {Q^ — Q^L{1 — k))~^ , while the Q'' observable is just —iNe~^Q''L. In Wegner's case 
the differences between and Q'' are even more daunting. And yet we will see later in this paper that Fyodorov's 
— formalism is able to reproduce and extend the SUSY results, including even the correct normalization factors. 

I submit that the desire for an obvious symmetry in the Lagrangian is misguided. Should supersymmetry mean 
symmetry between and Q^, or instead symmetry between these variables and the Q-^ — determinant? The 
Grassman variables were responsible for the determinants in the numerator; recall that when we integrated them the 
result was the product of the — Q'' determinant and {det{Q-f))^^^ . When calculating the two point correlator we 
will see that integration of and produces a matching determinant which almost cancels the determinant. 
The correct normalization factors come from cancellation between the two. 

Fyodorov's formalism possesses a hidden supersymmetry which is guaranteed by its exact derivation, protected by 
the Ward identities, and manifested in every detail of the Lagrangian. When calculating the density of states one 
finds that the saddle point is cx e"^, cx — e~"^; we see that the inverse relationship between the and 
observables preserves the Ward identities. The k in the determinant turns out to be crucial for extended systems, 
and its property A:|0) = is directly linked to the Wigner-Dyson statistics seen in the weak localization regime. Even 
the in [N — I^) \nx plays a role in obtaining correct overall signs. 

III. THE SADDLE POINT APPROXIMATION 

Analysis of this theory must begin with identification of the dominant parts of the Lagrangian and of the dominant 
values of and Q'^ . We will procede in two stages, first treating the eigenvalues x and only later treating the angles 
U and T. Concerning the eigenvalues, it is important to discern carefully which terms determine the saddle point 
values. A good first approximation neglects the determinants in equation 1181 in later sections we will take those 
determinants into account properly. Unlike previous authors, we do not assume that the saddle point is spatially 
uniform. The saddle point equations are 

= -7Vx„fc+A^(l-/''7V-i)a;;^i + 2ziV(l-/^Af-i)i/2sin0„fc, 

sin ^ 2-1(1 - i''N-'r'/\ul{r'Ef + iY^{{k/{\ - kyi^,u.,xiul)u.u 

sin0^ EE 2-\l - I^N-')-''\T-\r'& -zY,K.V2T.AT-')T.),, (24) 

Each component of x has two solutions controlled by a sign. We will encapsulate the sign signature of the saddle 
point in the diagonal matrix Lk = ±1. The solutions of the saddle point equation are x^yk 
LfcVl - /''iV-i cos {(j)yk) + iVl - I^N-^ sin (0^,^). The constraint x^L > imphes that = -^j- 

These equations describe the static equilibrium of + 1-^ identical particles moving in a quadratic potential divided 
by an infinite logarithmic wall. Each particle chooses to live close to one of the two minima of the potential, according 
to L. The position of each particle within the potential is determined by the external force 2sin0^/j. From a different 
perspective, these equations describe an energy band of width 4e \/l — I''N^^, and the position of 2sin0„fc within 
that band determines the phase 01,^. 

These saddle points dominate the model. All statistically significant configurations of and have x close 
to xq. Typical deviations can be estimated from the Lagrangian - at the saddle point its local quadratic part is 
— ^{x — xq)^, where rj^k = 2 cos (f).vkS.xp(—iLyk4'vk)- Imagining an added Langevin style of dynamics in the model, 
one sees that any deviations away from the saddle point values are subject to a restoring force equal to ~Nrj{x — Xq). 

Therefore fluctuations away from the saddle points will be of order {2N cos 4>) ' or smaller. The cos signals that 
the saddle point approximation breaks down when 2sin0 is close to the band edge. 

At first blush these solutions permit 2^^^ saddle points; i.e. a different L at each site. In reality (1) L docs not 
not fluctuate from site to site, and (2) only 1^ + 1 out of the 2 saddle points of L can be distinguished one from 
another. The latter point is the most accessible. Define P as the group of permutations of L which conserves Tr{L). 
P's only significance in the saddle point equations is to decide which eigenvalue a;£ is coupled to which force 2 sin (j)kv 
This choice is physically significant only at the global level, and then only if is not proportional to the identity. 
Therefore only a single global choice of P is required. This is a consequence of the theory's symmetry under rotations 
of [/. 

Going further, all members of P can be transformed into each other by finite rotations x^ —> UpxlU^. Therefore 
if one eventually performs a full (non-perturbative) integration over U then only a single member of P need be 
considered even at the global level. 
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What of the possibihty that Tr{L) might fluctuate from site to site? We will return to this question in section IVCl 
where we will establish that there is a single optimal value of L and that deviations from that value are penalized on 
a per-site basis by a free energy cost proportional to In k. 

The saddle point approximation breaks down when 2sin(/) is close to the band edge. The T variables in sm(f>'^j are 
unbounded, and therefore a complete integration of T always includes the band edge and beyond. Therefore we must 
either do the T integrations prior to the saddle point approximation, or else regulate T. This paper will rely heavily 
on the mechanism of spontaneous symmetry breaking, which regulates fluctuations in T, requiring them to be small. 

The global value of T is not regulated by SSB. Obviously this is not a problem if E'' = E'^, in which case Z contains 
information about self-correlations of single eigenfunctions. However one often wants to know about correlations 
of energy levels or eigenfunctions, and therefore chooses E^ 7^ in this case one needs some way to control the 
global T integral. There are two alternatives; the first is to perform the global integral prior to the saddle point 
approximation, and will be discussed in a little more depth in section fVI Bl The second alternative is to pin T. One 
pinning mechanism is to eliminate T {T = U = 1) by calculating just the density of states [I^ = 1^ — 1). Another 
pinning mechanism is to break T's continuous symmetry by choosing energy levels that are far apart compared to the 
level spacing A; i.e. \& — &\ ^ A. In the main we will depend on the latter pinning mechanism. 



A. Spatially Uniform Saddle Point and 1/g Corrections 

The saddle point solutions given in equation [M] arc not quite satisfactory because they are not spatially invariant: 
sin (jjk depends on U and T which fluctuate from site to site. Therefore we choose a spatially invariant saddle point 
and rely on perturbation theory to manage the difference between it and the correct saddle point. We extract the 
zero momentum component of U and T: Uy — > U^Uq, Ty TyT^, and then choose a saddle point which depends 
only on [/„ and Tq; sin(/){ = 2-i?-i(l - I^N-^)-^/^{uIeWq)^^, sh^ct)] = 2-ig-i(l - l''N-^)-^/^{T^^&To),j. The 
difference between this choice of sin (f) and the correct choice is handled during the integration around the saddle point: 



/ 



dx exp(-:^(a; - x^f + 2iN{\ - l''Ny/'^6{sm(j))) 



= Jdxexp{-^{x-xo-Xgy + ^xl),Xg = 2ir]-'^{l-l''Ny/'^6{sm(t>) (25) 



Therefore the perturbative prescription is simply to add Xg wherever x occurs (x Xq + Xg + x), and to add a term 
T X^ufc ''?fe^gt)fc ^'^ Lagrangian. These Xg corrections correspond to 1/g corrections, where g is the conductance. 
Although 1 /g corrections have been calculated in the supersymmetric literature^^, it seems likely that these particular 
corrections have been neglected because the standard SUSY procedure is to assume that the saddle point is spatially 
uniform despite fluctuations in U and T. This amounts to simply dropping Xg. 

The logarithm in the Lagrangian plays a crucial role in this perturbative approach. It generates cubic and higher 
vertices which are proportional to — 1)*^; the resulting perturbative corrections are well controlled as long as 

((^ — 1)^) < 1/N. Since x (x xo+Xg + {Nrj)"^^'^, the logarithmic diagrams are in control only when Xg/xo «C N~^/^. 

This is essentially a restriction that fluctuations in U and T, when multiplied by k, must be smaller than iV~^/^; i.e. 
a requirement of spontaneous symmetry breaking. Additionally, the size oi E ~ E is constrained. The perturbation 
theory breaks down and one must use the spatially- varying saddle points given in equation [M] if SSB is not observed 
or E — E is too big. 

With a little algebra the Lagrangian at the uniform saddle point can be simplified : 

- ^Y.^'ok + iN - Inxofe + 2iNVr\l - I^'N-^f' ^ xok sin</)ofe 

= -2-^NV{I^ + P') - A^V^^ sin^ (j}k + {N - ^ In(ife) + iLkicpk + cos</)fc sin</)fc) (26) 



B. Fluctuations in the Eigenvalues 

Having identified the saddle points, we begin a process of integrating out the eigenvalues x. Fluctuations away 
from the saddle point will be called Xy , giving the decomposition x„ = -I- Xgy + Xy . We decompose the Lagrangian 
into the part Co + Ci which does not depend on Xy vs. the part £ which does depend on these fluctuations. We also 
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approximate det(l — fc) w 1. 

Z = J J dUo dTo dU dT e^°+^^+^ 

Co = -iVy^sin2(/)fc+t(iV-/'')F^4((/)fc + cos0fesin0fe) 

fe k 



2 

vk 



~2 $Z ^^/^^ k))viV2Tr{Uy^XQUljJy^XQUl^) + ky^y^TriTy^x^T^^T^^x^T^^) 

ViV2 ViV2 

/ -n- 
dx e^' exp C^-^Vkxlk) det{a + H ^vdAii^ + ^iv + ^l) ^vdMi^o + ^gv + xl), 



V1V2 V1V2 
00 



vk 2=3 



Xok 



Lk\/l-I''N-^ e'^"'^\ T]k = 2cos0fc exp(-^Lfe,^fc) 
mct>{ = 2-ig-i(l-/*iV-i)-i/2([/t^/C/o).^, sin0^*=2-ig-i(l-/''7V-i)-i/2(ro-i#To),-,- 
= irii^~e-\ul{UlEfUy - E^)Uo + *e~^(/e/(l - k)U,UlU,,xlulU,)u « Af-V^ 

x^j = ir^-'r\To\T-'E^T, - E')To - ze ^ fc„„,T-iT„,a;^T-iT„),,- « 

V2 

N-1 -V 



^ „-2-iivy(/''+/-'')+2-iivye-^y^^''Tr(£;^£;-'),-i^JV/''-yjv/-'' 

^ ^ -zAre-i(T(4 + a;^, + x')T-'Uj, > ^ (27) 

rfJi- ^ dJt.dJt, dfi. dJ^ , ^ ' 

VJ Vljl V2J2 Vljl V2J2 

Several terms in these equations are controlled by 7V~^/^ and can therefore be neglected: these include the Xg +5^ 
term which multiplies the kinetic term within the determinant, the logarithmic vertices, and all instances of I^N~^ . 
Additional N~^^'^ terms occur within the Van der Monde determinants, which are multiples of many differences 

Xki — Xk2 ■ If a pair of signs L^^ and are the same then the corresponding factors of Xki — Xk2 nearly cancel; 
otherwise these factors are of order Lj.j2cos(6. Therefore we decompose the Van der Monde determinants into the 
multiple of two parts: a part Eq containing cosine factors selected by the condition Lk^ = — -Lfc^, and a second part 
E composed of the small factors selected by the condition Lk^ = ifcj • We neglect the instances of Xg + x occuring in 
Eo- 

e^' contains kinetic terms which are quadratic in x. These contribute to the Hessian, which is H ~ dPC/di:"^ = 
NVji^viV2^jij2 + where ^,;i«2jij2 = fcfif2(^u7^^f2)iu2(^ij^^^fi)j2ii- (The Q-f sector is similar.) The saddle point 
approximation generally requires that the entire Hessian should be included in the Gaussian kernel of the e^' integral. 
In our case ^ is controlled by the small parameter k, implying that the eigenvalue dynamics are almost local. Therefore 
we will treat (, perturbatively, leaving it out of the kernel and including it in e^' . The validity of this perturbative 
expansion could be analyzed easily if e^' were simply a Gaussian integral without the determinants and the logarithmic 
vertices: it would be a question of whether could be computed perturbatively. If fluctuations in T and U are small 
then ^ w fc <C 77 and perturbation theory is justified. The determinant and the logarithmic corrections complicate 
things; in this respect the perturbation theory's validity must be explored by calculating the first order corrections 
etc. and seeing whether they are small. 

In future expressions we will omit the integration over x, with the understanding that sooner or later this integration 
must be done, using Wick's theorem and {xyj^kiXv2k2) = SviV2^kik2^~^lk^ ■ 0-^ ^® were to include the full Hessian 
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then we would have (i^jfejiu^fc^) = A'' ^{r] + (,) ^.) However we already multiply the overall normalization constant 
by J]^(^)^/2. The simplified La grangian is 

= e^' dct{a + K) Eo{xo) E{xo + Xg, + x,) 

k V 



C. Correspondence to the SUSY Sigma Model 

We compare the saddle point action £i in equation [27] with the action of the SUSY sigma model^, 
^ J dx [—D'S/Q ■ VQ — 2tLuAQ]. Both Lagrangians contain kinetic terms which are quadratic in Q, and also mass 
terms which are linear in Q. The correspondence can be made precise by setting k — -jV • V, N'^^ = iru J dr, and 
neglecting the (1 — k)~^ factor multiplying Q^'s kinetics. This will allow us to reproduce SUSY calculations of the 
two point correlator and of anomalously localized states. 

However there is an important difference in the way that this paper's sigma model and the SUSY sigma model 
treat the energy band. The Lagrangian developed here incorporates band information explicitly. and Q'^ are 
proportional to p cx cos 0, and are small close to the band edge. Therefore the kinetics are more sensitive to the band 
edge than than the mass terms. 

In contrast, the SUSY matrix Q is conventionally normalized so that its eigenvalues are ±1; the SUSY model's 
only explicit band dependence is via the overall multiplier v. In other words, the SUSY model packages all energy 
band information in the diffusion constant D and the density of states v. These are understood as phenomenological 
constants. The present paper's Lagrangian can be construed as a derivation of D's band dependence: Du oc p^. 
However this interpretation mixes up two different kinds of physics which are probably best left separate: one is the 
Lagrangian's sensitivity to Q's fluctuations; i.e. H = d'^C/dQ'^, and the other is the fact that Q gets small close to 
the band edge. In any case, presentations of the SUSY sigma model in its final form typically do not provide explicit 
prescriptions for calculating how the kinetics and observables depend on E. 

To my best knowledge, this present article's sigma model is the first to include explicitly the kinetic term's correct 
band dependence. This explicit band information will be manifested in our results for the two point correlator, 
anomalously localized states, and perturbative corrections. 



D. Saddle Point Analysis of the Wegner Model 



Analysis of the Wegner model's saddle points obtains results quite similar to those of Disertori's model. We begin 
by calculating the spatially uniform saddle points. The following identities are true by construction: 

dZ dZ 
p{E) = -^-'L,Y,Im-^, p(E)=7T-'UY.^'^:7jr (29) 

Taking the derivative produces Re {Tyx\T~^) jj = Ljp{E), where p = p{E)TTeN'^V'^ . p scales with the ratio 
of the disorder strength to the band width; p oc {ek)~^ when the kinetics dominate and p (x 1 when the disorder 
dominates. This contrasts with the Disertori model's p ~ cos0 which is always of order 1. When evaluating the above 
equation's J-^ derivative we neglect l''N~^ and the Q-^ — coupling's dependence on - this term's contribution to 
p is surpressed by iV~^/^. We obtain Im{{iUx^U^ + ^k)~^)^^-- w Lip{E). In general the U and T averages prevent 

precise statements about Xq and Xq, but U and T can be pinned by calculating just the density of states {P = = 1) 
or by choosing energy levels that are far apart compared to the level spacing. Both mechanisms allow one to drop 
the averages over U and T and obtain Re {x^jj « Ljp{E), Im {{ix{ + ek)~^)^^ w Lip{E). In the weak localization 
regime (efc ^ 1) the latter equation allows us to connect p{E) to the density of states p^k of the kinetic operator ek; if 
Rex{ then Im {ix{ + ek)~^ = —'nV~^p^k{—I'mx{) sign{Re{x{)). Corrections to this equation are controlled by 
Rex{ /Ea, where E'a is the typical scale of variations in efc's density of states. (In a finite lattice x{ should be large 
compared to the level spacing in order to smooth the density of states.) We conclude that in the weak localization 
regime p{E) « TTV~^peki^~^E) and Im {x{) ~ —e^^E. When the disorder dominates (efc <C 1) x{ is never small and 
is probably best represented as a phase; one obtains Re {x{)~'-^ « —Lip{E). 
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This analysis can be made more precise by considering the equation governing spatially uniform saddle points . 
The saddle point equations are practically the same as Disertori's saddle point equations if the disorder is dominant; 
here we consider the weak localization regime efc ^ 1. We leave the — Q^ coupling and Van der Monde determinants 
out of the saddle point, just as we did with Disertori's model. 

= -Nxl^ + iNr^UlE^Uolu + tN{ixi^ + ek)^^ (30) 

When efc ^ 1 this saddle point equation is solved formally by x^.^ = iVi + p{^ri) sign{Re{xQ^)), where r.i — 

{UqE^ Uf))ii + Re{ek — ± . True solutions can be obtained numerically, and may yield non-perturbative 

information about the position of the mobility edge. We expect that the solutions will be qualitatively different 

depending on whether or not E is well within the energy band. The sign of Re{xQ^) is not fixed by the saddle point 

equations, signalling the existence of several saddle points distinguished by their sign signatures L. When E^ = E, 
f . . , . ... . . ' . . 

Xq 's imaginary part is proportional to the identity and its real part is proportional to L. This is the same structure 

seen in Disertori's model. 

Next we turn to saddle point analysis of Xq. This analysis is valid only when T is pinned; otherwise one must do 
the T integration prior to the saddle point analysis. 

= -Nx''^ + iNl-~\T^^E''To)jj + N{x''j)-^ + -j-^, e^"' = J dW exp{-i€NTr^j{SLk)) (31) 

We have seen that Cw must cancel the Q-^ logarithm's dependence on e, and that when T is pinned Re {x^)jj ~ 
Ljp(E); clearly the last two terms in the saddle point equation must combine to produce a resolvent like the one 
governing x^. Equation 1231 gives the leading order result for C\y in the e/c 3> 1 regime, which contains a logarithm 
— {N — I^/2)Tr In(Lefc). Corrections to £w must reproduce the resolvent by adding a self-energy to the logarithm's 
argument. As we have already outlined, calculation of these corrections is a rather intricate task and is outside the 
scope of this paper. When 1^ = 1 the first correction is proportional to N'^£XqJ2viv2 ^^i^2^viV2' ^ different 
formulation like the two-particle irreducible (2PI) effective actioni5'il2. niay be necessary. More general considerations 
indicate that when E'' — E, Xq^ can have only one of two values selected by the value of Lj. We expect that Xq 
will change qualitatively when E approaches the band edge. The biggest outstanding questions are whether Xq^s 
imaginary part is insensitive to Lj, and is the same as Xq's imaginary part. Both questions are answered affirmatively 
in Disertori's model, and also in Wegner's model when the disorder dominates. 

We conclude that the Wegner model and the Disertori model have the same saddle point structure: xq sa Lp{E) + 
is(L,E). The fermionic sector of s is proportional to the identity and is approximately equal to —e~^E. The same 
may be true for s taken as a whole. The main difference between the two models is p, which in Disertori's model is 
controlled by the simple semicircular energy band. In the Wegner model p is controlled by the spectrum of efc when 
efc 3> 1, and is the semicircular band when efc 1. 



1. The Hessian 

Fluctuations in the x-^ eigenvalues are controlled by the Hessian, the second derivative of the Lagrangian with 
respect to a;:^^: 

Eigenvalue fluctuations scale with the inverse square root of the Hessian; we see immediately that they are controlled 
by 7V~^/^, just as in Disertori's model. There the Hessian is NrjSy^y^ -\- Nk; it is almost local, which allowed us 
to decouple the eigenvalue integrals site by site. The Wegner Hessian is almost local both when efc 1 and when 
efc ^ 1. To establish this fact we must analyze the second term, which is the product of two Green's functions. It 
may evaluated numerically if accuracy is desired, but in the two limits efc ^ 1, efc ^ 1 analytical techniques may be 
used instead. In the efc <C 1 regime x{ is of order 1 and the two Green's functions may be expanded in powers of 
ek/x{; at leading order the resulting Hessian is the same as is found in the Disertori model. 

The opposite limit efc ^ 1 is less trivial. Direct examination of equation [321 is invalid because fc has eigenvalues 
that are arbitrarily small, including of course the zero mode fc|0) = 0. The standard analysis runs as follows. We 
switch to a momentum basis, in which the Hessian is 

Hy,y2isi,S2) ~NS,,,J{si - S2) - NS^.^Jisi - S2)V~^J2(''^( + efc(s3))-i(zx{ + ek{si - s^))-'^ (33) 

S3 
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The second term may be rewritten in terms of the density of states: 

- NV-^S,,,J{si ~ S2) J d<^p,k{<^){ix{ + <,)-^{{ix{ + ek{si ~ 53))"^), 

{{ix{ + ek{s, - S3))-i) ^ ^ - ek{h)){ix{ + , + (efc(si - S3) - eki^W (34) 

S3 

In order to make further progress analyticaUy one must restrict the momentum si to obey the condition ek{si — 
S3) — efc(s3) <^ Rexj Vs3 3efc(s3) = c;. This equality is certainly true when si = and may be true for a range 
of small momenta. The restriction means that our analysis will only give the long-distance behavior of the Hessian; 
shorter distances require numerical evaluation. Next we expand {ixj + ^ + {ek{si — S3) — eA:(s3)))~^ in powers of 
— (efc(si — S3) — <;)/{ix{ + ^). In the efc » 1 regime Rexj/EA <C 1, equation [Ml mav be converted to a contour 
integral, and one obtains 

H^iV2{si,S2) = -NV^^6i^iJ{si - S2){V + 27ri^^^^^^^ - Tn-^{p^k{'i){ek{si - S3) - <;)))^^_„/ (35) 

The first two terms are mass terms, while the last is kinetic. Rex{ « p « nV'^pek oc (efc)~^, so the first term 
is roughly (efc)^ larger than the last two terms, and the Hessian is nearly local. This same reasoning provides the 
standard derivation of the low energy Lagrangian governing Goldstone bosons both in our models and in the SUSY 
sigma model. In these cases the mass terms are either absent or exactly cancelled by another term in the Lagrangian. 
Therefore the Goldstone bosons are massless and their low- momentum Lagrangian is proportional to ek oc Z3V • V. 

The obvious difficulty with this reasoning is that as ek becomes larger our ek{si — S3) — ek{s3) ^ Rex{ condition 
becomes ever more restrictive, and for large enough ek becomes the same as the simple restriction si — 0. One must 
ask at what point does the restriction make the entire analysis meaningless; how small can we make the low-momentum 
Lagrangian's range of validity? It is not clear whether moving to the continuum really answers this question. Probably 
equation 1351 should be interpreted as a reliable analysis of orders of magnitude at small momenta, but not as a recipe 
for calculating either the Hessian or the Lagrangian. If either is required one should evaluate the Hessian numerically. 

Turning to x'', its Hessian is proportional to N, and therefore its fluctuations are controlled by iV~^/^. Analysis of 
the Hessian's locality requires computation of the second derivative of Cw- We saw that Cw^s first derivative must 
produce a resolvent, so it is likely that its second derivative is a product of two Green's functions, the same as in the 
x-^ sector. 



2. Corrections to the Sigma Model Approximation 

Here we analyze what happens to the saddle point when U and T are not spatially uniform. As with Disertori's 
model, we do not allow the saddle point to fiuctuate but instead introduce a perturbation theory in Xg = H~^d, where 
d is the difference between the spatially uniform saddle point equation and the true saddle point equation: 

Hxl = iNi-\U^o{UlE^U^ - Ef)Ua)^, + iN{U^ {tUx^ U'' + ek)-^U)^„, - ^iV(^4^ + e^)™ (36) 

The N's on the right side cancel the N in the Hessian. If fiuctuations in U are small then the last two terms may be 
expanded in a Taylor series; one obtains a form analogous to the Disertori model's Xg, which is 

4v^ = ^V~'e-\U^oiUlE^U, - Ef)Uo + »e^(fc/(l - k)UMUv,xlulU,)u (37) 

Similar results can be obtained for the bosonic correction Xg . The vertices of the Xg perturbation theory correspond 
to third and higher derivatives of the Lagrangian's logarithms and of Cw- These vertices are proportional to N, so 
the Xg perturbation theory is in control only if Xg <C N~^^^. Any non-locality of Wegner's Hessian should not change 
this requirement significantly. Because of the complete correspondence between Wegner's Xg perturbation theory and 
Disertori's perturbation theory, all subsequent convergence analysis applies equally to both models. In particular, we 
will find that the sigma model approximation can not be controlled except in the spontaneously broken phase where 
fluctuations in U and T are small. 
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3. The Sigma Model Lagrangian 

We write Wagner's sigma model Lagrangian in terms of a part £i which controls spatial fluctuations of U, T and a 
part £o which depends only on Uq,Tq: 

Ci = iNe-^ Tr{xiul{UlEfU„ - Ef)Uo) + iNr^ ^ Tr{xlT^\T~^E''T, - E'')To) 

V V 

+ {N - l'')Try,{ln{Ux{,U^ - lek ~ irKjf) -lii{xl - lek)) + CwiJ'') - Cwi-J'' = 0,T = 1) 
£o = ^^Y.''ik+^NVr'TriU„xlu^^Ef)+tNVr'Trinx'„T^'E') 

vk 

+ (TV - /'')Tr„(ln(4 - + (iV - Z'') ^1 ^^^vj + '^w{j'' = 0, T = 1) 

S,,,,_L = Cls,,WlWy,s,,C,,L, Clsla,L = T,xlT-\ e^'^' ^ J dW cxpi~iNTr,j{SLiek + r\P))) (38) 

The saddle point equation ensures that £o's leading dependence on E - E is just iNVe-^Tr{UoxlU^Ef) + 
iNVr^TriTox^oTo^E''). 

Unlike the SUSY sigma model, the present sigma model is not fully described by the Lagrangian. Much of the 
interesting physics lies in the — determinant, and also one must take into account the Van der Monde determi- 
nants, as described in equation!^ We have glossed over the — Q'' determinant's dependence on W, which causes 
a coupling between the determinant and the sigma model Lagrangian. 

In the efc ^ 1 regime this sigma model is capable of producing accurate predictions incorporating full information 
about the energy band. The logarithm controlling C/'s kinetics can be expanded in a Taylor series whose coefficients 
may be determined either numerically or via analytic approximations. Cw^s Taylor series coefficients are equivalent 
to certain moments of the dW integral, and these also may be calculated numerically or approximately. An important 
outstanding problem is whether Ituxq is the same as Imx^. If so, the kinetic terms and the E terms in Ci are 
respectively quadratic and linear in p, the same as in Disertori's model. If not, in subsequent formulas one must 
substitute for p. Later we will calculate observables in Disertori's model; the reader should bear in mind that the 
Wegner Lagrangian is very similar and produces analogous results. 



IV. THE - COUPLING IN THE WEAK LOCALIZATION REGIME 

We now narrow our focus to weak localization phenomena; i.e. the regime in which the kinetics dominate the 
disorder. More precisely, in Disertori's model the weak localization regime is characterized by the condition^S that 
7V-1/2 <^ ko,k <^ 1. Exh = efco is the Thouless energy, the smallest non-zero eigenvalue of the kinetic operator ek. 
We will use the small parameter (TV^^^fco)"^ to control the theory. 

We will also assume spontaneous symmetry breaking of the continuous global symmetries of U and T, implying 
that the system is in the delocalized phase. In I? = {1,2} dimensions the Mermin- Wagner theoreroii prohibits SSB in 
very large volumes; however in finite volumes one can still see an effective SSBiS,. We will outline a way of calculating 
corrections to the SSB assumption in powers of the inverse conductance, preparing the way for rigorous proofs about 
observables in the weak localization regime. 

Q-^ and are coupled by a determinant det(a -I- k), where 

Ct = ((a^Q -|- a^guj + ^vi)ii ~ (-^O '^gvi ^vi)ji)^iii2^jij2^viV2j ^ — (2^0 ^" ■^gvi ^ ^tii )fcDiD2'^fi-^t)2 ^viUl^ (39) 

Our assumption of spontaneous symmetry breaking implies that fluctuations in U and T are small; as a result k's 
spectrum is controlled by k. In particular k has a zero eigenvalue fc|0) = 0; therefore A = a + n has near-zero 
eigenvalues that are particularly sensitive to fluctuations in x, U, and T. In order to capture this physics we will 
project out the zero-momentum sector of A = a + n and treat it separately from the other momenta. We define two 
projection operators: Pq selects out the zero- momentum modes, and P+ = 1 — Pq- The following formulas describe 
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the decoupling: 

det(A) = det{P+AP+) dct{PoAPo - PqA{P+AP+)-'^APo), 

PoAPo = (4 + 4 + iO-(xg + 4 + i'') + i^-i4^(fcTr-i;7[/^),,,, 

V1V2 

v'/\PoAP+u.,., - (4,^ + i4 - 4 - - + - 4 - 

+ x'o Y^ikTT'^UU^,,, - V-^xl {kTT-^UU%,,, (40) 

V3 V3V4, 

The new x notation prescribes spatial averaging y-i The kTT-^UW terms are exactly zero when U and T are 
constant; they are proportional to the fluctuations in U and T, the same as Xg. If fc is long-ranged then the sums over 
V are regulated by factors of V~^^'^ caused by the absence of long range correlations in U and T. 



A. The P+ Sector 



In the following paragraphs we will analyze the spectrum of P^AP-^ when Xg = x — 0, starting with a's spectrum 
and later taking k into account, a is a local operator; its eigenvalues are just Sxoij = — x^y Xq^ and nearly 
cancel each other if i and j satisfy Li = Lj ; in the opposite case they add to each other. The small eigenvalues of a are 
quite small, of order 6xq cx {Ei — Ej)/ [i-qe) if Ek doesn't depend much on k. The large eigenvalues are approximately 
Li2 cos{E) . We separate the two sectors of a with projection operators Psmaii + Pi = 1; the projection operator for 
the small eigenvalues is Psmaii = Ps = ^i^t^^titi^hh^^ivi- 

The — coupling includes an additional term k which mixes the eigenstatcs of a. We estimate the effect of mix- 
ing on the large eigenvalues by applying the identity a\^) = — F+;A(F+sAP+s)^^y4P+;)|V-'). The 
term is of order ^^2 cos(i?) -l-O(fc), while the last term in the identity is of order x (Pj^sAPj^s)^^ . We have assumed 
that fluctuations in U and T are small; therefore P+s^-P+s max{k, 6x0); in consequence k makes only perturbative 
changes in the large eigenvalues, with the changes controlled by the small parameter k. Similarly, we analyze the 
effect of the mixing on the small eigenvalues with the identity a\ip) = {P+gAP+s — P+sA{P+iAP+i)~'-^AP+s)\'ip)- The 
P+sAP+s term is of order max{k, (Sa;o): while the last term is of order fc^, so its effects on the spectrum are controlled 
by the small parameter k. We conclude that P+AP+ has a sector well described by P+j whose eigenvalues are all of 
order Li2cos{E) and another sector well described by P^s whose eigenvalues are all of order fc. Corrections to this 
picture of the spectrum are controlled by the small parameter fc, and are very small compared to the saddle point 
action (equation [57]) which contains a factor of N, so we will neglect them. 

Because the smallest eigenvalue in P+AP+ is of order ko, the effects on its spectrum of Xg, x, and fluctuations in U 
and T are perturbative, controlled by (iV^/^fco)^^ ^ 1- Expanding the P+AP+ determinant perturbatively produces 
terms like exp(Tr(fc~"'^(a;g -f- x))). These terms are very small compared to the saddle point action, so we will neglect 
them, approximating det(P+j4P+) as det(P+s(a;Q — a;Q -f- a;Qfc)P+s) det(P+i(a;Q — XQ)P+i). 



B. The Zero Momentum Sector 



We turn to the determinant which controls the zero momentum sector, det(PoAPo — PoA{P+AP+)~^APo). We 
introduce projection operators Pq;, Pqs satisfying Pq = Pqi + Pqs and then analyze the determinant in terms of the 
two sectors, g"^ signifies the kTT~^UW terms. For simplicity we assume that these terms are smaller than N^^^^. 
The following matrices give orders of magnitude; the upper left entries correspond to the Pq; sector while the lower 
right entries correspond to the Pqs sector. 



1 9-' 
g-^ 6x0+ N-^/^ + g^^ 



PqAPo « 
PoAPo - PoA{P+AP+)-^APo(x 



, PoAP+« 



9 



-1 



, (P+AP+)-i« 



1 1 

1 fc-i 



N-^/^+g 

g-^ - g-^N-^/^k-^ 6x0 + iV-1/2 + g-i _ AT-ifc-i - g-^N'^/^k'^ 



(41) 



We have dropped terms whose relative magnitude is N~^/'^ or fc or smaller, and also some terms of order g^^k^^ . In 
the second line the matrix elements coupling Pqs with Pq; are less than iV^^/^, while the Pqs — Pos matrix elements 
are of order N~^/'^ + 6xq. Therefore the coupling between the two sectors is less than A^~^/^ and can be neglected. 
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producing det{PQi{xQ—XQ)Poi) det{PosAPos — PosA{P+AP+)~^ APqs) . The contribution to the Pqs determinant from 
mixing with the P+ sector is controhed by {N^/^ko)~^ , which is small in the weak localization regime. We neglect 
this contribution, arriving at a final form for the — coupling in the weak localization regime: 

det(a + K) = det(P+a(4-4 + 4fc)^+«) det(P,(4 - 4)^0 det(Pos^Pos) 
Po{a + K)Po = (4+x^ + i^)-(a;[;+4 + i'') + y-i4^(/cTr-iC/C/t)„^„, (42) 

ViV2 

We have derived this form under the implicit assumption that we will not be calculating correlations between Q-^ 
and Q''. The — determinant is the only source of such correlations in this theory, and should be treated more 
gently when they are being computed. 

C. Decoupling of IJ and T 

The mathematically controlled reasoning leading to equation 321 is extremely interesting because it sketches a 
proof that \J and T decouple almost completely in the SSB phase: the only exception is the zero-mode coupling 
det(Pos^Pos)- Consider the case of 1^ = 1-^ = 2, which is appropriate for computing two point correlators. If U and 
T were decoupled completely, the U sector would be just the classical Heisenberg model, while the T sector would be 
a hyperbolic sigma model introduced by Spencer and Zirnbauer. The hyperbolic model exhibits SSB unconditionally 
in D > 2; there is no phase transitioni^i. The Heisenberg model in D > 2 has been proven rigorously to exhibit 
spontaneous symmetry breakinglS,. It is also believed that these models display an effective SSB even in _D = {1,2} 
when the system is sufficiently small. 

The Q-^ — zero-mode coupling contains only a few factors multiplying the path integral, while the saddle point 
action (equation [27]) is proportional to the volume NV. Roughly speaking, the Q-^ — coupling is a factor of {NV)~^ 
smaller than the action, and is best treated as a prefactor. This implies that the zero-mode coupling does not change 
the SSB behavior of the underlying Heisenberg and hyperbolic models, and that the SSB assumption is self-consistent 
in Z? > 2 dimensions. 

In order to prove SSB in Z? > 2 dimensions one would have to prove bounds on the spectra of P+s^P+s and 
P+iAP+i which are stringent enough to show that P+AP+ is bounded below O{ko). Such bounds may be available 
from probabilistic arguments based the fact that the theory penalizes configurations in which the determinant has 
one or more small eigenvalues. 

D. The Wegner Model 

Wegner's — coupling is a bit different than Disertori's: 

deti-iek + Qf - A\S)), = " i^k),,,,,^,,{Qf - iek)-X,,,Jv,v2nj2Lj„ 

ii.vi 

(43) 

Ignoring fluctuations in W, at the saddle point S is spatially constant and the determinant simplifies to det(— lefc -f 
Xq — a;o). This should be compared to the Discrtori model's det(a;Qfc + ~ Xq). We see again the same possibility of 
cancellation vs addition between Xq and Xq, which distinguishes the Pg sector from the Pi sector. When the disorder 
is small (efc ^ 1) the saddle point solutions are roughly RexQ « p oc (efc)~^, so the Pi sector's eigenvalues are not 
much larger than the P, eigenvalues. The only real distinction between the two sectors is that the P+; eigenvalues 
have a real component of order 2p while the real component of the P+s eigenvalues is much smaller and perhaps even 
zero. The real component corresponds to a mass, and it is remarkable the mass in Wegner's Q-^ — determinant 
are rather small compared to the kinetic energy (2/5 + lek), unlike the Discrtori model where the P; sector's masses 
completely dominate the kinetic energy. 

The rough equality of the P+j and P+s eigenvalues simplifies our argument that none of the P+ sector's eigenvalues 
is much smaller than kg, which is the basis our factorization of the Pq sector of the determinant from the P+ sector. 
Within the Pq sector Pqs can still be factorized from Pq; , so that a controlled analysis of the — Q'' coupling in 
the efc ^ 1 regime should obtain a form similar to Disertori's coupling given in I42[ and eventually arrive at identical 
predictions for observables at leading order. 
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V. FLUCTUATIONS IN AND 

When the global symmetries in and Q'^ are spontaneously broken, fluctuations in these fields are small and 
weakly interacting. The next step is to integrate these fluctuations. To simplify the arithmetic, we specialize this 
paper's mathematical formulas to calculate two point correlators. We set = = 2, so that Q-^ and Q'^ are both 
2x2 matrices. The dominant saddle points, i.e. the ones which minimize the number of near-zero factors in the 
determinants, are Lf = (0-3)4^, = {'^3)jj for the Retarded- Advanced correlator, and L/, — {1,1,— 1,-1} for the 
Advanced- Advanced correlator. (0-3 is the Pauli matrix.) At the — {1,1,-1,-1} saddle point all of the U and 
T terms disappear from the saddle point Lagrangian (equation [57]) leaving only X^ufc '?fe^gt)fc; ^^i^ saddle point 
describes almost local dynamics. Therefore in 13 > dimensions the Advanced-Advanced correlator is at first order 
purely local, the same as the zero-dimensional result. We will focus on the Advanced- Retarded correlator. 

We adopt parameterizations of T and U which are tailored for perturbation theory: 



U = 



tyf - zf ^l-(y/)2_(^/)2 



(44) 



The inverses of T and U can be obtained by inverting the sign of both y and z. y and z vary from ~oo to 00, while 
y^ and z^ vary within the unit circle defined by (y-^)^ -I- {z^)'^ — 1. The integration measures are dy^dz^ and dy^dz'' . 

We parameterize Uq with \l,Ol such that Aq = 1 — 2{ylY — 2{zlY , tan^Q = z^/y^. Similarly we parameterize 
To with Xl = 1 + 2(2/^)2 + 2{z^)^, tan0^ = z^/y^ The Jacobian for Uq^Tq in these coordinates is 2-^dXld\ldel^del 
The limits of the Ag integration are from 1 to 00, while the Ag integration is from —1 to 1. In these coordinates the 
following relations are very useful: 

iU^UiasUlUi)!! = ~{UlUia3UlU2)22 



A2(l - 2y( - 2zt) + 2Jl-XlJl-yl- zfiyi cos 6*2 -I- zi sine's) 



{2yl + 2zi - \)(2yl + 2x\ - 1) + 4^1 -y\-zl^\- yl - z|(yiy2 + ^1^2) 
1 - 2(yi - y2f - 2(zi - Z2f + 4(j/2 + z\)[yl + zl) - 2{y\ + z\ + yl + zl){yiy2 + ziZa) 
2)22 



A2(2yf + 2zf + 1) - 2./A2 - Ia/I + yl + z^yi cose^a + 21 sine's) 



= {2yl + 2zl + l){2yl + 2x1 + 1) - ^^l + ylVzi^l + y(+4{yiy2 + 21^2) (45) 
« 1 + 2(yi - 2/2)' + 2(zi - Z2f + 4(y2 + zl){yl + z|) - 2(y2 + zl + t/^ + z|)(yiy2 + ^1^2) 

We define E = E + {u!/2)a3, uj = Ei — £'2. The saddle point equations are 2esin(/)/ = (JJgiJ-'j/o)^ = + 
iw^A''(cr3),,, 2esin(/)5 = {T-^&Tq)^, ^E+ ^uj''X''{a3)jj. Their solutions are 

- E 

xo = CT3(e*'"^'^-f z(2?77fe)~iAcj) + 0(w2) = a3p + s, sin^= — 
p = cos^ -I- 4- 0(ti>2), s = tsinff) ~ ^t&ncf) + 0{uj'^); 



Co = iNV^Lki<j)k + cos (t)k sin (t)k) - NV^ sin 



= -NV2-^r^{{Ef f + {E'' f) + iNVr^{io''X^ cos^^ + uf XUos^f) ~ NV2-^r^{{uj^X^ f + {u^ X^ f) + 0{uj^) 
= -NV2-^r^{{Eff + {E'^f) + iNVr^{Lo''X^ + X^ pf) + NV2-'^r^({uo^X^f + {Jxff) + 0{Lo'^) (46) 

Under the assumption that loe^^ ^ 1, we drop the quadratic term in the last equation, set Tr{Ef Ef) - {E^f - 
{E^Y ~ 0, set p = cos(/), and set (2^)^(1]^ ?7fc)"^/^ » 1- We also use the Stirling approximation [N - 1)\{N - 2)! w 
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27rAf2iV-2g-2iv simplify; 



N-l 

l=N-2 k 

X exp {-2NV + 2-^NVe-'^Tr{EfEf) - 2-^NVe-'^{{Eff + {E^f) + NVie-^{uj^X^p + uj^X^p)) 

« 2-4iv2^7r-2^y2y(2p)2V' {iNVr^u^'' X"" p + uj^ p)) e'^' d.et{a + k), 

det(a + K) = (2/5)2^ dct(P+,,(i(2?r/fe)"i(A^w'' - A'^cj'') +x|;cr3fc)P+s) det(PosCr3(a + k)^^o,s) 

Pos(J3{a + k)Po« = «(2f77fc)-i(A^a;^ - A^'o;'') + ^3(4 + 1^ - ^3(0;^ + ^*') + V-^xl<Jz Y.i'^TT-^UU^U^^ (47) 

In calculations of the two point correlator Pg merges the i and j indices: P, = SiijiSiii2^jij2^viV2- In some places we 
have neglected the difference between E^ and E^, but this can be restored easily if needed. 



A. Integration of the Fluctuations in U and T 

We have already assumed spontaneous symmetry breaking, implying that fluctuations in U and T are small. We 
now integrate the small fluctuations, reducing the model to zcro-dimensional integrals. For the moment we ignore e^' 's 
dependence on these variables, and will return to this issue later. Earlier we moved the zero-momentum component 
of U,T into Uo,To; therefore J2v — ^^'^ terms in the Lagrangian which are linear in y,z come to 

exactly zero. The equations arc further simplified by noting that the part of xq which is proportional to the identity 
does not play any role in U and T's dynamics. To second order in y, z the integrand and action are 

£1 = -2iNcofX^pe-' Y.((yvf + + '^^NuJ'>X'>p~e-' J^{{ytf + {ztf) + 0{iNujrWl - X^{y\ z^)) (48) 



_ (^/(l - k)Uv2{yLyL + 4A2) - -yp'T. kv,v2{ylX2 + <4) + 0{Nky\ Nkz'') 

V1V2 V1V2 

±- ^ -^iVe-l(^^ + cSa3(A''(l + 2y^y^ + 2zU) - 2^ (A'')^ - l(y„ cos^g + z, sin0g)(l + 0{y^ z')))) 
(P d d ^ 



dJy^ 



"75b TJb ' lib lib ' — P + {x^gvl — X^gv2 + ^tl ~ ^v2)/2, Sy — S + (x^^i + x''gy2 + 5^1 + ^v2)/2 

aJvijx""Jv2n "'•''vijl ""^V2j2 

We neglect the (1 — k)~^ controlling Q-^'s kinetics and perform the Gaussian integration, which generates a mul- 
tiplicative constant equal to {AN p^ / 'k)'^~'^'^ , plus two determinants. The c subscript on the expectation value in the 
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following equation denotes the fully connected part of the expectation value. 



X 



det"^(fcP+ + iP+ujf \f {2pi)-^) det"^(fcP+ - iP+uj''\^{2piy 



= (det(a + n) e^'+^ '"'^s-'')c/,t , 

det(a + «;) = dct(F+,(«(2?77fe)-i(A-''w^ - A^'lj'') + x[;a3fc)F+.) det(PosfT3(a + '«)^0s) 

Po«CT3(a + k)Pos = «(2e~77fe)"'(A/cj^' - X''lu') + aaC^ + F - - I'') + 0{g-^) 

— — (^/(-^ ~ ^))viV2'^'^{^Vl^Vi Uv2^V2 ) ~^ "o" ^ v ^l^i t'2 ^'^(^t'l '^V2^V2'^V2 ) 



i;ii;2 i'i'U2 



{yiyj = (8iVp2)-i^„^|n(_(2^g)-i^'>A'')|«2>(l + 0(.9-2)), 

«y4) = (8iV/52)-i(i;i|n((2p6~)-ia;/A0|«2)(l + O(5-')), 0(7) ^ (P+fc + P+*7)-i (49) 

n (conventionally called the diffuson propagatopS^i^,) describes propagation of only non-zero- momentum states. We 
will verify that p — NVpn^^e^^ is the density of states. Equation [49] is valid only if SSB occurs; in particular if SSB 
does not occur then we expect {yviyv2) 1. Moreover this equation is only the leading order result, and is subject to 
corrections: 



The kinetic term Suiua ^viV2Tr{T^^a'iT~^^T^^azT^^) contains every even power of y'', z''. [/'s kinetics are 
similar. These terms generate perturbation theory vertices with even numbers of legs; every vertex is proportional 
to Np^k. The vertices, which we will call l/gp corrections, modify the free energy density and change the value 
of correlators hke {yviUvi}- (The conductance g is defined^S a.s g = Eth/^, where A is the level spacing. In 
Disertori's model g = pko.) 

• The observable -jj^ contains every odd power of y^, z^, causing more 1/gp corrections. 

• The mass terms generate vertices with every odd power of y, z, and all are proportional to Npute'^^. Like the 
1/gp vertices, these mass vertices modify both the free energy and observables. 

• Xg occurs in many places, both in the Lagrangian and in observables. It contains a kinetic term which is 
proportional to rj^^pk and contains every even power of y, z, and also a mass term which is proportional to 
r]^^e~^uj and contains every odd power of y, z. 

• Cj contributes additional terms with every even power of y, z, multiplied by Nk{xg + x)"^ . 

First we consider the 1/gp corrections. All of the vertices are proportional to Np^k, while the bare propagator 
2/112) °^ {vi\{N p'^ k)~^\v2) ■ Assuming that in perturbation theory diagrams each vertex with its associated Np^k 
roughly cancels a single propagator {yviyv^)^ diagrams with I loops are controlled by {yviyv^) ~^ oc {Np^k)^~^ . This 
inverse assumption breaks down in 13 = {1, 2} where SSB does not occur— - the diagonal elements of the propagator 
diverge in the V ^ 00 limit. However, in finite volumes one still has [N fP'Y~^ scaling, so for large enough N p^ one 
obtains a controlled perturbative expansion. 

As long as w(2pe)~^ ^ k, the mass vertices are much smaller than the 1/gp vertices. However once uje^^ is of the 
same order as the conductance g the mass vertices become more important. At the same time the low-momentum 
behavior of the propagator changes to cx {AN p uj / e)~^ . In this regime the Feynman diagrams are controlled by powers 
of [AN pu}/e)^^ instead of 1/gp- 

We can now check the validity of the SSB assumption and the expansion in powers oiy, z. The leading 1 /gp correc- 
tion to the free energy density is proportional to Np'^V^^ ^vrv2 kviV2{XviXy2)'^ oc {Np^k)^^. As long as this diagram 
is small compared to 1, we may conclude that the perturbation theory is justified, and that the SSB assumption is 
correct. In fact we have already seen that in D > 2 dimensions SSB does occur in the weak localization regime. Even 
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when D = {1,2} one stih sees an effective symmetry breaking in small enough volumesiS; the free energy density 
correction is a good gauge of how small is small enough. The in the correction is a sign that close to the band 
edge the 1/gp perturbation theory breaks down and fluctuations become large. 

Turning to Xg^ we have already seen that Xg/xa <^ N~^^'^ is a necessary condition for control of the logarithmic 
vertices; if this condition is violated then the spatially uniform saddle point for the eigenvalues is invalid, and one 
must use a spatially fluctuating saddle point. We can now estimate the magnitude of Xgj xq - its kinetic part is of 
order kl^y^^y.-o^ . If SSB is observed then this is oc {Np^)~^ <^ 7V~^/^; however if SSB does not occur then the kinetic 
part is of order k which may be much larger than N^^^^ . Estimating the magnitude of Xg/x^s mass term requires 
a little more thought, since it contains odd powers of y, z and therefore must always be paired with another Xg or 
with a mass term. The simplest estimate is x^/xq cx {2pe)^'^ uj^ {yy-^y^^)] in the SSB regime one has the condition 
we-i < 24p3. 

Lastly, Ci is proportional to k times all even powers of y, z. As long as one is not calculating correlations of x, its 
effects are smaller than the other kinetic terms, and may be neglected. 

B. Correspondence to the SUSY Sigma Model 

The supersymmetric sigma model has been used to make detailed calculations of corrections to leading order results; 
see for example Blantner and Mirlin's work on correlations of eigenfunctions [2l| . These corrections were in powers 
of 1/5, and correspond to our 1/gp expansion. The extra p in this paper's perturbation theory signals the fact that 
we have explicitly included band information while the SUSY model does not. 

In the present model both £/ and the Xg terms arise from fluctuations in the eigenvalues; probably SUSY papers 
have not calculated them until now. It seems likely that they have confined themselves to the present model's 1/gp 
vertices and mass vertices. 

C. Fluctuations of the Saddle Point Signature 

Our saddle point analysis indicated the possibility that each site would separately choose its own saddle point 
signature L, leading to a sort of Potts model. Further analysis of this scenario seems hopeless, since one would have 
to integrate the U fluctuations in the presence of a disordered background corresponding to the various site-wise 
configurations of L. Here we will establish that there is a single optimal value of L and that deviations from that 
value are penalized on a per-site basis by a free energy cost proportional to Infc. 

Consider the two determinants in the denominator of equation I49i which were obtained by integrating out the U, T 
fluctuations. They were obtained for the Retarded- Advanced saddle point L = {1, —1, 1, —1}, which is equivalent to 
the { — 1,1,1,-1} saddle point. Now consider the other Retarded- Advanced saddle point L = {1,1,1,-1}; at this 
saddle point the U fluctuations are local and do not produce a determinant in the denominator. The difference means 
that the former saddle point is exponentially favored over the latter saddle point by [k + i{2p)~^uj\i~^)~^'^^ . 

Now consider site- wise fluctuations away from the favored L = {1,-1,1,-1} saddle point. If a single site switches 
to the disfavored saddle point then the entire path integral is multiphed by k + i{2p)^^oj\e^^ <C 1. This per-site 
penalty controls the fluctuations and avoids the Potts model scenario. 

This argument can be generalized easily to other saddle points. When calculating the Advanced- Advanced correlator 
the favored saddle point is i = {1,1,-1,-1}. I have checked all saddle points required for three and four point 
correlators {P = =4). The general rule seems to be that the favored saddle point minimizes \TrL\, and that the 
per-site free-energy cost is roughly i(TrL)^ Infc. These results were obtained for the Disertori model; the Wegner 
model requires more careful analysis before determining whether and how its saddle point signatures are regulated. 
This is likely the only major remaining conceptual challenge concerning the Wegner model's weak localization regime. 

VI. OBSERVABLES 

The supersymmetric sigma model has been used to calculate many different observables. Here we calculate just 
a few, to check on agreement between the two models, and to finish our analysis of how to maintain mathematical 
control. We will drop Xg and 1/gp corrections. Blantner and Mirlin^-'^ found that while 1/g corrections are important 
for eigenfunction correlations, they are do not enter at lowest order into the level correlator R2. 
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A. Large u) approximation 



Equation |49] has the following form: 

/+1 roc I* 

d\f j j dOfdO^ e"P^'''^'+^'^'^ Zx{\^ A\6f ,e^) (50) 

Z\ was calculated using a saddle point approximation which fails when lo^X^ is large enough to move the saddle point 
to the band edge. Our solution to this problem is to "pin" A'' by making uj^ large. This does not actually restrict the 
d/^ integration; the real logic is as follows. We assume that the Z\ is much larger inside the band than elsewhere, and 
that the proper way to control the A'' divergence is by losing information about the band edge. Therefore we perform 
a Taylor series expansion of Z\ in powers of A** — 1 and then integrate. When TipLo^ 3> 1 the exponential oscillates 
very quickly, so that the d\^ integral is dominated by the A = 1 limit of integration. In this scenario the band edge 
really is unimportant, and we obtain a systematic expansion of Z controlled by powers of (Trpcj'')^^. 

One may easily verify that this procedure is equivalent to changing coordinates to A** = 1 + 2(yQ)^ + 2(zq)^, 
expanding Z\ in powers of y and z, and then performing the resulting Gaussian integral. These are exactly the 
same coordinates which we used to integrate T's spatial fluctuations; from a certain viewpoint we have simply done 
a Gaussian integration of all of the T degrees of freedom together, including Tq- Based on this viewpoint, Andreev 
and Altshuler— claimed that they had applied the saddle point approximation to all of the T degrees of freedom 
including Tq. This explanation is a bit misleading, since properly speaking the saddle point approximation is a Taylor 
series expansion of the exponent, while in the Tg integration the exponent is treated exactly at leading order in 1/ gp. 
The real Tq physics is in the Taylor series expansion of Z\. 

Andreev and Altshuler— used this Taylor series expansion to perform both the f7o and Tq integrals. When 
calculating two point correlators (J** = 1^ — 2) there are two saddle points corresponding to the two permutations of 
e(,E2', in the angular coordinates which we have chosen these permutations are represented by Ay = ±1. Altshuler 
and Shklovskii^^ computed the A-^ = — 1 contribution nine years earlier than the Andreev- Altshuler publications. 
They avoided generating functions or sigma models and instead resummed perturbative expansions of the Green's 
functions, finding a smooth oj~'^ behavior. Later Andreev and Altshuler used the supersymmetric sigma model to 
calculate R2 completely, including both the smooth part and an oscillatory part coming from the A/ = 1 saddle 
point. Here we will reproduce and extend their results by applying the same saddle point approximation to the present 
sigma model. 

The new coordinates are A** = 1 + 2{y^)^ + 2{z^f, = ±(1 - 2{y^)'^ - 2{z^)'^), where the ± specifies the saddle 
point. The yo,zo integrations should be done at the same time as the yv,Zv integrations. We revise equation [49] to 
match the altered procedure; 

Z = det-\k±iLof{2pe)-'^) det-\k~iLo\2pe)-'^)e'^^P'~''''^''^'^^^ 

X det{P+s{t{2iriky^{\f'ojf - X'^lj'') + xla3k)P+s) det(i(2gr;fc)"^(A^w^' - X^u^) + a^ii^ - i^)) 

^ ^ -^N^\sl - {cl{l + 2{{ylf + {zlf)))) 



-N'r\si + + 2{{yt,r + «)')))(4 - 4(1 + nvir + 



d d 7'' 

+ ^N'e-'tltlMfiy'j' + i<n<?)c (51) 
(«> = {8Np')-\v,\ik~^{2p~e)-'u;'')-'\v2), (y^y^O = (8iVp2)-i(t;i|(fc ± z(2pg)- V)->2), 

This can be simplified a bit, dropping x(y^) terms because they are 1/gp smaller than the x terms, and factoring 
out the zero modes. 

Z ^ A:tp^{±ujfJ')-^d&X-^(kP^ ± iP^J{2plY^) det-\feP+ - iP^J'{2pl)-^)e'^^P^''^'^'^'^'^^^ 
X det(iP+(2g77i)-i(±cj-^ - w^) + P+(/5 + s)fc„,„J det(iP+(2e?72)~^(±tj-^ - tj'') + P+(/5 - s)fc„,„J 
X {^(2l^^)-^{±uJ^ - J") ^'x{~x\) {i{2lT]2)-^{±uj^ - J") -x{^x\) 

^ ^ -^N^\xo2 + il - 2p{{ylf + {zlf)) 



d P 

a2 

^_^^2~~2i I ~h I oA/„,fc„,fc I J>\\ ('^__ I zi> OA/„,f>„,fc I J> J>\ 

jb J jb 



-N'r\xoi + il^ + 2p{ylyl + z^zD) {x^^ + - 2p{ylyl + z^zD) 



d T° dl° 

+ 47V'?-'p'(y.^y^,+«)c (52) 
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Next we do the x integrations, set to-^ = lu'', and make the uji~^ ^ 1 approximations 771772 = 4/5^, Xq^Xqj = 
{p+s){—p+s) « —1. The density of states can be determined using p = tt^^ J^v Im{dZ / dJy^j^2), and at leading order 
is given hy p — NV pii^^e^^ . We have speciaUzed the formahsm to calculate the two point correlator (/•'' = 2), 
and in this case p seems to be dressed by corrections in powers of 1/gp. These \/gp corrections must completely 
cancel, because the density may be calculated in the — ~1 case where there are no angular variables. However 
the density may be dressed by Xg corrections to the saddle point and by terms proportional to k from the determinant. 
The second derivative is 

T. 9 9 / „ N dct^(2pefcP+) 



det(2/5efcP+ + iujP+) Aet{2pekP+ - iujP+) 
+ {{2pek - iLu)-l^^f + N^r^ - i2Nr^uj-^p{l - i{2pl)-^uj)-^ {v\{k - i{2p~e)-^uj)-^\v) 
+ V~'^{Tr{2plkP+-iujP+)-^f (53) 

Next we subtract the second derivative at the Advanced- Advanced saddle point. As we mentioned earlier, this 
saddle point does not involve any kinetics at highest order; it produces the zero-dimensional random matrix result, 
— iV^e^^e^^*"^. To obtain the two point correlator i?2(w), we take the real part, sum over vi and W2, and divide by 

2^2p2. 



i?2(w)-l = {27:^p^)-'^ReTr{{2pek-iuj)-^) 
+ (27rV^)"^cj"2cos(27rptj) 



det^(2p?fcP+) 



det(2pefcP+ + iu)P+) det(2pefcP+ - iljP+) 
+ {2t:'^ p^r^ Re {Tr{2plkP+~iLjP+)-^f, p ^ cos4> + 0{uj) (54) 

To make contact with the SUSY result, we choose the kinetic operator to be k = — -jV^. Except for the last term, 
our result for i?2 is the same as that obtained already by Andreev, Altshuler, and Shklovskii^ i^^i^^i^^ , but with added 
information about the energy band: we have a new p multiplying the kinetics. The kinetics are more sensitive to the 
band edge than the mass terms, and as a consequence i?2's oscillatory component is amplified near the band edge. 
To my best knowledge, this present article is the first derivation of the two point correlator which explicitly includes 
this energy band information. 

The oscillatory term characteristic of Wigner-Dyson level repulsion would not exist if the eigenvalues obeyed the 
Gaussian statistics one would naively expect of massive modes in a sigma model; in that case the x\^^ and £21)2 
the observables would integrate to zero. Instead the Q-^ — determinant's factors of x^ — x^ distort the probability 
distribution of the eigenvalues, causing a repulsive force between x^ and x^ which pushes x^ toward positive values and 
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toward negative values. The origin of Wigner-Dyson level repulsion is in eigenvalue repulsion'^^. This non-trivial 
eigenvalue dynamics is hidden in the SUSY approach where eigenvalues are integrated quite early, but it is still there, 
working through the Grassman variables. 

The oscillatory term comes from the A = -fl saddle point, and has a A = — 1 sister which contributes to the 
last term in P2- This sister term is of the same order as the other terms, and is controlled by the trace of j/'', z^'s 
propagator. It was not reported by Andreev and Altshuler, who may have believed that a factorizable term can not 
contribute to R2{uj) — 1. Despite its appearances it is produced in a non-factorizable way: part of it comes from the 
(^1111^21)2-^1^2) expectation value, and can not be decomposed into {x\jj^x\x2) x (^2^2 ^1^2) =0x0. Verification of 
whether the extra term is a bona fide improvement of Andreev and Altshuler's results would require evaluation of 1/(7 
corrections which might be able to cancel it. The extra term's structure can not be reproduced by l/g corrections 
to the prefactors or the determinant, which leaves perturbative corrections from the Lagrangian as the only available 
cancellation mechanism. Evaluation of such corrections is outside the scope of this paper. 



1. The Wegner Model 



The ratio of determinants controlling i?2's oscillatory term is particularly interesting because it is sensitive to all 
of the theory's degrees of freedom. The determinants in the numerator originate in the P+s sector of the — 
determinant, while the determinants in the denominator come from integration of the the U,T degrees of freedom. 
The eigenvalue degrees of freedom also contribute to the denominator, but in Disertori's model these are exactly 
cancelled by the P+/ sector of the — coupling, so the net contribution of the massive modes is a factor of 1. 
This agrees with the SUSY treatment of the sigma model approximation, in which the massive modes integrate to 
exactly 1. Our partial analysis of the Wegner model leaves open the possibility that massive modes may contribute 
non-trivially to that model, especially if the w-induced mass of the P+; sector of the — coupling could become 
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larger than its small natural mass 2/5 <C 1. In this case one would see a knee in the oscillations' magnitude around 
the point where oj becomes a significant contribution to the eigenvalue mass. Understanding this issue would require 
a careful solution of the Wegner saddle point equations and evaluation of Wegner's Hessian and determinant. 

B. Small uj approximation 

When to is not much larger than the level spacing the band edge is important, or in other words the complete 
non-compact Tq manifold must taken into account. Therefore one must perform the Tg integration prior to the saddle 
point approximation. The only term in the Lagrangian f equation llSp which depends on Tq is 

iNi-^ Tr{T,ToxlT^^T-\& - .1'')) 

V V 

- iNr^^X'^X-^Y.^^' - ^Jv){<i - 42)^1 + 2/^2/^ + 44 [coselvl + sin^gz^) (55) 

V 

Integration of Tq at an early stage would result in an infinite-range interaction with all powers of y, z, causing problems 
for further analysis. Therefore Kravtsov and Mirlin^^ began by integrating the spatial fluctuations, producing a zero- 
dimensional model with special terms in the Lagrangian reflecting the fluctuations. Afterwards they integrated the 
spatially uniform component of the SUSY matrix Q over the complete domain of integration. Kravtsov and Mirlin's 
results are complementary to Andreev and Altshuler's results; corrections to the former are well controlled when 
Loe~^ <C 2/5fco, while corrections to the latter are well controlled when [NV p)~^ <C Loe~^ . Together they provide a 
complete description of the two point correlator for all values of uj which are well within the band edge. Kravtsov 
and Mirlin's technique has been reused for calculating many other observables, including eigenfunction correlations, 
multi-fractal statistics of eigenfunctions, and \/g corrections^^. Here we will reproduce their two point correlator by 
integrating [/q and Tq over the whole domain of integration. 

Kravtsov and Mirlin's approach requires doing the saddle point approximation last, and therefore is difficult to 
reconcile with the controlled approach we have adopted in this paper. The saddle point approximation supplied 
information which allowed us to control the — determinant and other issues, and was a cornerstone of our 
integration of the spatial fluctuations. This is not a difficulty in the conventional SUSY derivation, which is not as 
attentive to control issues. For us it is a difhculty. It might be possible to simply assume the saddle point results, 
perform the integrations, and then justify the saddle point assumptions a posteriori. More likely the correct approach 
is to start with the Tq integration and then find some rigorous way of controlling the global interaction, perhaps by 
treating it separately from the rest of the theory. Here we will not pursue such possibilities but instead imitate the 
SUSY approach, foregoing mathematical control. 

The main purpose of this calculation is to demonstrate that the — models can reproduce the standard SUSY 
two point correlator at small lo. In particular we are neither looking for additional terms like the one we found in the 
previous section, nor searching for any other improvements. Such efforts belong to a more considered and rigorous 
treatment. Kravtsov and Mirlin's result depends only on the four point correlator {yt-^vl^}- Since we only want to 
reproduce Kravtsov and Mirlin's result, we will drop all higher order correlators and all occurences of the on-site two 
point correlator (y^). 

We start with Disertori's model, as given in equation[TSl and develop it as follows. Using the Stirling approximation, 

the normalization constant becomes 7 = 7V^^2^^^7r^'*^e ai^"^' '-^^•'e^^^. It is necessary to delay the integration of 
spatially uniform values of x^ and x'' until the last moment, but the spatial fluctuations in these variables can be 
integrated immediately, multiplying the path integral by Wki'^^)^^'^ ~ (""/^p)^^"^- The angular variables are 
accompanied by a 2^'^ Jacobian. We assume that x^,x^ are close to the Retarded- Advanced saddle point Lk — 
{1, 1, —1, —1}, and therefore set x{ — x^ = x\ — X2 — 2/5 everywhere except in the spatially uniform part of the mass 
term. We use the weak localization form of the Q-^ — determinant given in equation [42l We neglect corrections 
to the saddle point approximation, i.e. Xg and £/, set & — & .lo^ — iJ' ^ neglect l''N~^, and make several we ^ 1 
approximations. The resulting path integral, prior to the next steps, is probably equivalent to the SUSY sigma model. 

Turning to the — Q^ determinant, we assume that feg 3> — |, |a;2 ~2;2|} and we expand the the P+s determi- 
nant in powers of x^ —x'', keeping the second order contribution and dropping all other terms. (These steps require that 
we-i < 2/5fco.) The determinant becomes -{2pf^V^ det^{P+k){x(-x\){x^-xl^) exp{-2-^Tr{P+k-^) Y.^ixi)~'^ (4 - 
x^)^). Lastly we expand U and T in powers of y, z, and neglect all cubic and higher powers; these correspond to 1/gp 
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corrections and may be easily restored if needed. The result of these steps is 
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Next we integrate the spatial fluctuations y,z. Because we have omitted all cubic and higher terms this Gaussian 
integral may be done exactly, but we follow Kravtsov and Mirlin and treat the lo terms as perturbations, assuming 
again that bje~^ <^ 2pko. In any case these terms must be treated perturbatively when we do the A integrals. The 
integration multiplies the path integral by tt^'^"^) (2/5)~^^^~^)7V~^(^~^Met^^(P+A:), and correlators are determined 
by (2/1,1 2/^2) = p'^)~^ {i'i\{P+k)~^\v2) and Wick's theorem. To determine the effective Lagrangian induced by the 
integration one calculates fully connected diagrams with no dangling legs. As announced earlier we keep only diagrams 

proportion., to the fourth „Le„. ia.vt^ The diag™ w'h L fae.o. o, is ex.c% 1 because 
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We first outline the correct procedure for doing the angular integrals, and then take instead a shortcut. One should 
first turn the terms which are quadratic in A into derivatives with respect to sources, then do the remaining trivial 
integral, and then take the derivatives, as follows: 



e^bd'^/df 



/=0 



(za + Z)-ie 



1 ia-\-l 



1^0 



(58) 



Applying this technique to equation [571 one produces two factors in the denominator: {l^ + mpiu) — 6Jy))^^{l^ + 
nTpu!)~^, p = NVpe^^TT^^ . The Vs and J's in these factors dress the observables, multiplying the number of terms 
which should be evaluated. We simply neglect these terms, which are not necessary to reproduce Kravtsov and Mirlin's 
result. Applied to equation 1581 neglecting the sources in the denominator produces the final result ia~^ exp(ia + b). 
This procedure is justified only when a = npuj 3> 1, but one may verify that in the zero-dimensional case the real 
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parts of the neglected terms completely cancel, so there is no effect on the D = final result. 

Z = 7VV22V-2p4g2^-2^2ivygi^T.(££) y^^ ( dx\dx\dx{dx{ [x{ - x\){x{ - x\) e''+^^" 

^ Jx'>L>0 

C^fj = -2-'p-'r' + (c. - W,^)(c. - S,Pj){v^\{P+k)-'\v2f - 2-'Tr{P+k-')Y,ix',)-'{x{ - x\f 

^ = -^Y.'^l + NVY,^nxu--iNVr\j\x\ + ?^x\) 

k k 

+ iNVr'^{{E + uj/2)x\ + {E~ oj/2)xl) + iNVr^{{E ± uj/2)x{ + {E T t^/2)4) (59) 

Next we take derivatives with respect to the sources. Corrections to the first derivative may be dropped because 
in the 1^=1^ = 1 calculation of the same quantity there are no angular variables and therefore no corrections; 
—iNi^^Xj. Taking the second derivative produces a term which is proportional to lu and can be neglected 

because it is purely imaginary, and another term proportional to which we neglect because it contains the fourth 
power of (ui|(P+fc)~^|u2). The remaining terms are 

2-^e-^p-^{vi\iP+k)-^\v2f -N^l-^x\xl (60) 



(7 7'' (7 7'' 

Lastly we perform the saddle point approximation. We neglect the effective Lagrangian's contribution to the saddle 
point equations, which is supressed by N~^. The saddle point solutions are x\ = exp(i(/)-|_), = — exp(— «(/)_), a;{ = 
exp{t(l)±),xl = -exp(-i0zp), where sin0± = (i? ± w/2)/(2e). The saddle point action is -2NV ~ NVE'^e^'^ + t{l ± 
l)puj. We shift x^ and x'' to match the saddle points; x{ — x^ — > — xJ + i(cr3)ii(l=Fl)w(2e?7i)~^. After performing this 
shift we neglect the remaining x^ — x^ occuring in C^jf and simplify; in the + term C^fj — —uj'^{2pe)~'^Tr{P^k~'^), 
while in the — term it is zero. 



Z = N^V^2^7:-^p'^i^Lo-^y2Te''-^^^'^'"^ [ dx\dxldx{dxl (4 ' 4 + T l)iUj{2eT]i)-^) 
X (4 -x\-{lT l)^w(2e772)"^) exp(£e// - 2-^NV^r^kxl) 

k 



(7 7" d T' 



2-'r'p- {vi\iP+kr \v2r - N'r'{eMi(^) + exp(-i0) + 4) (61) 



2V2 



The remaining Gaussian integrals multiply the result by 7r^iV~^/5~^. The terms which are proportional to uj are purely 
imaginary and may neglected. We expand exp(£e/^) « 1 + Ceff, drop a term proportional to the fourth power of 
{vi\{P+k)~^\v2), and simphfy. 

p = TT-^Y^ Ira — - = TT-^ ^ Im{iNr^ exp(-z0)) = NVpn^^r^ (62) 

V V 

= {2lp)-^v^\{P+k)-^\v2f - V-^u-^ + N^r^ - V-'e''P'^i2ip)-'Tr{P+k-') + V-'io~'e'''"^ 



d 1° (7 7' 



2V2 

To obtain the two point correlator we subtract the Advanced- Advanced second derivative — iV^e^^e"^"^, take the real 
part, sum over vi, V2, and divide by 27r^p^: 

R2{iu)-l = (2Tr^p^)-\cos{2puj) ~ 1){lu-'^ ~ {2ip)-'^Tr{P+k-'^)) (63) 

We choose k — -jV^ to make contact with Kravtsov and Mirlin's resulij^i, which is identical to ours except for the 
extra band information contained in p. 



C. Anomalously Localized States 



Until now we have considered only spatially uniform saddle points. Other saddle points which are not spatially 
uniform may also be important - for instance instantons, skyrmions, etc. Muzykantskii and Khmelnitskii studied 
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non-uniform saddle points in the supersymmetric sigma model and found that they control the conductance in the 
long time limit^^. Their non-uniform saddle points correspond to statistically rare eigenfunctions whose weight is 
concentrated in a small volume - these are called anomalously localized state o^^i^^'^^ . In the weak localization regime 
ALS control the asymptotic tails of the probability distributions of many observables; see Ref. for a review. Here 
we will show that the ~ Lagrangian produces the same ALS already seen in the SUSY sigma model. We will 
derive the same equations controlling the ALS which Muzykantskii and Khmelnitskii used as the starting point for 
their calculations. 

We begin by fixing the saddle point and parameterizing [/, T, = UxlW , and Q''L — TxqT~^ as 



U = 



isinip^e 



T = 



cos{2Tpf) 



-zsin(2?/''^)e 
-cos(2V'^) 



coshV'*' sinh'0''e 
sinh'0''e'^ coshi/j'' 
f 



s + p 



cosh(2V''') 
sinh(2i/''')e* 



- sinh(2V''')e-' 
- cosh(2?/;'') 



(64) 



The inverses of U and T can be obtained by adding tt to 9^ , 9^. varies from 1 to oo, ij^^ varies from to 7r/2, and 
9f 9'' vary from to 2tt. The Jacobian is i sin(2-0'^) sinh(2-0^), and this must be multiplied by another 1/2 because 
the Q-f parameterization covers its domain twice. 

We neglect contributions to the Lagrangian (equation [27]) from the Jacobian, from fluctuations in the eigenvalues, 
and from the (1 — k)^^ factor controlling Q^'s kinetics. We have seen that in the weak localization regime the — 
coupling reduces to a zero-momentum coupling, so we neglect its effects on the ALS. Lastly, we set E-^ — & . The 
remaining Lagrangian is 

£ = zA^e-i^rr(Q{^)+zA^e~-i^rr(Q^L£;)-^^fc,,,,Tr(Q4Q4) + ^^fc,,„,Tr(Q^^LQ,t,^) (65) 

V V ViV2 ViV2 

Following Efetov, we choose fc = f V • V, NJ2y='^''I dr. The La grangian becomes 



+ 



dr 
■KvD 



^^p^{2S/^^) ■ (2VV''') - iTTi^pujr^ cosh(2V/) + ^^p^{2\/'i{jf) ■ (2VV'^) - ini^pcor^ cos(2V'^) 



drsiT?{2i}jf)\/9f ■ \I9f + smh^{2il/')\/9^ ■ \/9^ 



(66) 



Clearly —C is minimized if 9^ and 9^ are constant, in which case the last line of equation 1661 is exactly zero. —C 
should be compared to Muzykantskii and Khmelnitskii's equation 7: 



A=^ J dr{[D{V9)^ -2iwcosh6'] + [^(Ve'i)^ + 2zw cos^i] 



} 



(67) 



Their 9, 9i, and uj should be identified with this paper's 2-0'', 2ip 



f 



71", 



and 



respectively. Except for the 



band information contained in p, the two equations are identical. Furthermore Muzykantskii and Khmelnitskii's 
Str{AQ) = 4(cos6li -cosh6') is equal to this paper's j{Tr{Qla3) +Tr{Qlc73)) = 27/5(cos(2V'-^) - cosh(2V''')), where 7 
is a positive multiplicative constant. They impose the boundary conditions 9 = 9i — at the interface with an ideal 
conductor, and it seems reasonable to conclude that the corresponding boundary condition 9^ = 0, 9^ = tt/2 should 
be used in this paper's model. With these correspondences in place, every result in Muzykantskii and Khmelnitskii's 
article on ALS^^ follows from this paper's model. 



VII. ASSUMPTIONS AND MATHEMATICAL CONTROL 



Having completed the process of calculating observables, we now review our assumptions and the small parameters 
which we used to maintain mathematical control. The really weighty control issues were centered on two tasks: 
deriving a sigma model via controlled integration of the eigenvalues, and controlling the Q-^ — determinant. 

Everything rests on the assumption that a good saddle point can be obtained while naively neglecting the — 
determinant and Van der Monde determinants. (One can perform a saddle point analysis including the determinants 
by moving their logarithm into the Lagrangian and then performing a Taylor series expansion. While the results 
are qualitatively correct, N~^l'^ no longer regulates the eigenvalue integration and one loses mathematical control. 
Mathematical control requires that the determinants be treated separately from the saddle point approximation.) 
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Inclusion of the determinants in the integral must not cause changes larger than N~^^'^ in the mean value of x'',x^ . 
Our later analysis justified this assumption in the weak localization regime. Outside of this regime two mechanisms 
protect it: (1) multiplies the saddle point action but does not occur in the determinants. (2) The determinants 
always increase the free energy; they create repulsive not attractive forces. The true ground state is likely to be one 
that minimizes the free energy associated with the determinants, and therefore should have some clear connection 
with the naive ground state. 

Our controlled derivation of the sigma model was based on six small parameters: 

• {Nri)~^ <^ 1. Consequently the sigma model breaks down when E is close to the band edge. As a further 
consequence we are required to control the non-compact T integrals. This is only possible in the SSB regime. 

• (ttwp)"^ <C 1. This constraint controls the non-compact integration of T's spatially uniform component. If it is 
violated the band edge becomes important and one must integrate Tq prior to the saddle point approximation. 
This second approach works only when uj{2pETh)~^ 1. 

• Spontaneous symmetry breaking; (?/t,i 1/1,2) 1/.9P ^ 1? validity of perturbation theory in y,z. This mechanism 

controls the non-compact integration of T's spatial fluctuations. 

• The diffusive regime, fc <C 1. The condition 1 — k > was essential to the Hubbard- Stratonovich transformation 
of the fermionic sector. Additionally the small parameter k allowed us to approximate the eigenvalue fluctuations 
as being local, and to control deviations from locality with perturbation theory. These issues were packaged 
in £7. We eventually argued that £7 does not contribute at leading order in N unless one is calculating 
eigenvalue correlations. While eigenvalue locality may not be absolutely essential to mathematical control, it 
circumvents a host of complications concerning the integration of eigenvalue fluctuations within the determinants 
and observables. 

• w(2pe)~^ ^1, fc ^ 1. Our naive saddle point analysis predicted several saddle points at each site corresponding 
to different values of the saddle point signature Tr{L). In order to avoid a Potts model scenario we had to argue 
that all sites share the same global saddle point. Our argument rested heavily on the assumption that lo and k 
are small compared to the band width, in which case the determinant resulting from integrating C/'s fluctuations 
selects a single optimal saddle point and imposes a per-sitc energy penalty for deviations from that saddle point. 
If a; or fc were of the same order as the band width then this argument would collapse and we would have a 
Potts model. 

• Xg/xo <C iV~^/^. The previous small parameters completely justify a saddle point integration, but the saddle 
point is not spatially uniform; it depends on the configuration of the angular variables U and T . Wanting 
to maintain contact with the SUSY sigma model and to avoid rather daunting mathematical complications, 
we made a further approximation and chose a uniform saddle point. The difference between this saddle point 
and the true saddle point can be managed with perturbation theory only if the two are reasonably close to 
each other. Xg, which quantifies the fluctuating forces caused by U and T, must be small. This requirement 
requires spontaneous symmetry breaking and constrains u). If these requirements are not fulfilled then one loses 
mathematical control and must revert to the non-uniform saddle point, obtaining a different and much more 
complicated sigma model. 

Having derived the sigma model, we turned to the task of controlling the Disertori model's Q-^ — Q^ determinant. The 
small k and ui parameters are the cornerstones of this analysis, allowing us to identify two sectors in the determinant 
and establish that the two don't mix significantly. There was also an auxiliary assumption of spontaneous symmetry 
breaking. The small parameter N~^, in conjunction with our choice not to calculate — correlations, allowed us 
to neglect the determinant's contributions to the action. With this foundation in place, we introduced a new small 
parameter: 

• Weak Localization Regime: {N'^/'^ ETh/^)~^ = (N^^'^ko)'^ 1. In this regime the determinant is dominated 
by the kinetic operator k; it simplifies to a zero- momentum coupling responsible for Wigner-Dyson level spacing 
statistics, plus a normalization constant. 

It was also convenient to take {gpN^^'^)~^ <^ 1, but this assumption is not necessary for control of the — Q'' 
determinant. 

In the weak localization regime the — model is dominated by the physics of U and T. Our SSB assumption 
gave birth to a machinery for calculating observables. We saw that in D > 2 dimensions the SSB assumption is self- 
consistent in the weak localization regime, since the — model is (excepting prefactors from the zero-momentum 
sector of the determinant) a product of two theories which have each been rigorously proved to exhibit SSB. In 
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D = {1,2} dimensions SSB is expected to still occur in small enough volumes, and we have provided a formula for 
estimating the numerical accuracy of the SSB assumption in any particular system. 
During calculation of the two level correlator we made additional assumptions: 

• We dropped all instances of Xg and 1/gp. These corrections can be computed easily enough, as demonstrated 
by analogous SUSY calculations of 1/g corrections. 

• We again exploited u!e~^ <C 1, A; <C 1. These are conveniences introduced for the purpose of analytical treatment, 
and could be easily avoided at the cost of complicated formulas or numerical integration. 

VIII. OUTLOOK 

The SUSY sigma model is widely understood as an effective field theory describing completely the long-wavelength 
physics of disordered systems. While this model has had its greatest successes in the weak localization regime 
and in the exactly solvable cases of = {0, 1} dimensions, it also gives a picture of the more strongly disordered 
conducting regime and of the Anderson transition. Its picture describes the weak localization regime in terms of 
spontaneous symmetry breaking and small fluctuations of the order parameter. In one or two dimensions the Mermin- 
Wagner theoreniii guarantees that fluctuations will become large in the V ^ oo limit, corresponding to Anderson 
localization. In _D > 2 dimensions the Anderson transition is expected to be analogous to the SSB phase transition that 
occurs in other D > 2 systems; at low temperature/disorder SSB occurs, while at high temperature/disorder SSB is 
frustrated. Within the framework of the SUSY sigma model, this intuition suggests that under renormalization group 
transformations (i.e. increasing the volume) the model will evolve from the weak localization regime into one of two 
stable fixed points: a localized fixed point and a delocalized fixed point. The disorder strength determines the flow's 
starting point, while the SUSY group structure and the system geometry control the flow. The SSB/renormalization 
picture is complicated by the expectation that the order parameter is a function'^'', and by debate about whether 
scaling is controlled by a single parameter— i^i^i^^. 

This present work gives a different picture of disordered systems. First of all, the sigma model approximation 
breaks down when SSB is frustrated, including at the band edge. More important are our results about the weak 
localization regime in D > 2 dimensions. We have seen that in this regime renormalization group transformations 
do not under any conditions evolve into a localized fixed point, which is especially striking because the same theory 
reproduces and extends the SUSY weak localization results, including anomalously localized states. This conflicts with 
the conventional wisdom that ALS in a given dimension D are precursors of localization; i.e. signals of a competing 
localized fixed point existing in the same dimension D. We conclude that while the SSB/renormalization picture is 
able to give an account of localization in _D = {1, 2} dimensions and near the band edge in higher dimensions, it does 
not give the full story about localization in the bulk of the spectrum. 

The most novel feature of the — models is the spectral determinant coupling with . We have seen 
that the weak localization regime is characterized by the condition that this determinant be dominated by the kinetic 
operator k. fc's smallest eigenvalues scale as |s|^ oc V~'^^^] as the system size is increased their control of the 
determinant will cease. In the V ^ oo limit the low end of the determinant's spectrum will be controlled by Q, not 
k. This suggests a new mechanism of localization in the bulk of the spectrum: the determinant's small eigenvalues 
may amplify Q's long-wavelength fluctuations. 
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Here we develop the unitary ensemble characterized by Hermitian matrices. This work can easily be generalized to the 
orthogonal ensemble^ and probably also the symplectic ensemble. 

In addition to the condition Nx^ek 3> 1, W's spontaneous symmetry breaking may also require that Q*" exhibits SSB. 
The SUSY literature generally has a looser definition of the weak localization regime. We have chosen this strict definition 
in order to preserve mathematical rigor. 

This is likely the reason why replica calculations of the two point correlator are unable to obtain the oscillatory component 
unless they take into account replica symmetry breaking. 



